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Preface 


The  analysis  of  the  dyneunic  response  of  pnevimatic 
transmission  lines  has  been  the  subject  of  an  ongoing 
series  of  studies  at  the  Air  Force  Institute  of  Technology 
School  of  Engineering.  This  thesis  is  the  latest  in  this 
series.  The  purpose  of  the  study  has  been  to  examine 
specific  areas  in  an  attempt  to  answer  questions  raised 
in  the  earlier  studies.  The  emphasis  has  been  on  the  con¬ 
tinuation  of  work  conducted  by  Captain  P.  G.  Wilkins  in 
1984. 

I  would  like  to  take  this  opportunity  to  thank 
several  individuals  without  whose  assistance  this  study 
would  not  have  been  possible.  My  thesis  advisor. 

Dr.  M.  E.  Franke,  provided  invaluable  guidance  and  advice. 
Dr.  W.  Elrod  and  Captain  W.  Cox  also  contributed  welcome 
assistance  in  specific  areas.  Captain  Phil  Wilkins,  whose 
earlier  thesis  provided  a  basis  and  direction  for  this 
study,  gladly  answered  my  questions  and  offered  insights 
into  the  subject.  I  would  also  like  to  take  this  oppor¬ 
tunity  to  express  a  special  thank  you  to  my  wife  Suzanne 
for  her  encouragement  and  confidence  which  has  made  my 
studies  at  AFIT  so  much  easier. 


—  Stephen  J.  McKenna 
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Abstract 

The  theoretical  predictions  of  the  dynamic 
response  of  circular  cross-section  pneumatic  transmission 
lines  were  compared  with  experimental  results.  Tests  were 
conducted  in  blocked  line  situations  as  well  as  in  both 
laminar  and  turbulent  mean  flow  conditions.  The  lines 
used  in  the  study  were  24  in.  in  length  with  inside 
diameters  of  0.041,  0.119,  and  0.195  in.  Tests  were  con¬ 
ducted  over  the  frequency  range  20-1000  Hz  and  over  a  range 
of  mean  flow  Reynolds  numbers  from  0-16000. 

Theoretical  solutions  were  obtained  from  computer 
simulations  using  a  laminar  flow  model  and  two  turbulent 
flow  models.  These  models  were  developed  from  the  Navier- 
Stokes  equation  and  include  components  to  account  for 
viscosity  and  radial  heat  transfer  in  the  fluid. 

An  examination  of  the  turbulent  models  was  con¬ 
ducted  based  on  a  nondimensional  frequency  and  attenua¬ 
tion.  This  led  to  modification  of  one  of  the  turbulent 
models  to  include  a  frequency  dependent  term  in  the  con¬ 
ductance  . 

The  effect  of  flow  conditions  upstream  of  the  test 
line  and  the  effect  of  orifice  design  on  the  dynamic 
response  were  studied. 


Experimental  results  showed  good  agreement  with 
results  of  earlier  studies.  The  use  of  a  bypass  device 
upstream  of  the  test  line  was  found  to  have  no  effect  on 
line  response. 

The  reconfiguration  of  the  orifice  devices  to  more 
closely  match  theoretical  models  produced  no  noticeable 
change  in  the  results. 


DYNAMIC  RESPONSE  OF  FLUID  TRANSMISSION  LINES 


I.  Introduction 


Background 

The  behavior  of  periodic  pressure  signals  in  fluid 
lines  has  been  the  topic  of  many  studies  in  the  past.  The 
results  of  these  studies  are  important  in  applications 
such  as  pneumatic  and  liquid  control  systems  in  aircraft, 
jet  engine  test  facilities  and  measuring  systems.  Both 
analytical  and  experimental  studies  have  been  conducted  in 
an  effort  to  develop  accurate  models  for  the  prediction  of 
the  dynamic  response  of  fluid  filled  lines. 

As  a  signal  travels  along  a  fluid  line,  the  magni¬ 
tude  and  phase  angle  of  the  signal  are  modified,  depending 
on  such  parameters  as  the  line  length,  diameter,  signal 
frequency  and  fluid  characteristics.  The  ratio  of  the 
downstream  magnitude  to  the  upstream  magnitude  is  referred 
to  as  the  gain  and  the  change  in  the  phase  angle  is 
referred  to  as  the  phase  shift.  The  ability  to  accurately 
predict  the  gain  and  phase  shift  is  of  major  importance  in 
many  system  designs. 

Studies  have  shown  that  the  model  developed  by 


Nichols  (Ref  1) ,  and  modified  by  Krishnaiyer  and  Lechner 


(Ref  2)  accurately  predicts  the  response  of  signals  in 
lines  with  blocked  and  laminar  mean  flow  conditions. 
Attempts  to  extend  the  use  of  this  model  to  turbulent 
mean  flow  conditions  have  shown  good  agreement  at  high 
frequency  but  result  in  lower  attenuation  than  observed 
experimentally  at  low  frequencies.  Studies  by  Brown, 
Margolis  and  Shaw  (Ref  3) ,  Moore  (Ref  4) ,  Briski  (Ref  5) 
and  Wilkins  (Ref  6)  have  concentrated  on  the  development 
of  a  model  which  can  accurately  predict  the  response  under 
turbulent  flow  conditions.  Briski  (Ref  5)  and  Wilkins 
(Ref  6)  worked  with  the  model  developed  by  Moore  (Ref  4) 
to  produce  modified  models  which  provide  accurate  and 
relatively  simple  means  of  predicting  line  response  under 
turbulent  mean  flow  ouiiuiLions. 

Objectives 

The  following  objectives  were  established  to  study 
the  dynamic  response  of  pneumatic  transmission  lines. 

1.  Review  current  models  to  ensure  proper  applica¬ 
tion  of  theory  to  experimental  conditions. 

2.  Verify  experimental  procedures  and  extend 
results  obtained  in  previous  studies. 

3.  Investigate  the  effect  of  conditions  upstream 


of  the  test  line  on  line  response. 


4.  Revise  existing  computer  program  to  incorporate 
revised  models  and  improve  documentation. 

5.  Investigate  effect  of  improved  orifice  termina¬ 
tion  devices  on  correlation  of  theoretical  and  experi¬ 


mental  results. 


i 


II .  Theory 

Electrical  Transmission 
Line  Analogy 

The  theory  of  electrical  transmission  lines  is 
well  developed  and  may  be  found  in  many  references.  This 
theory  may  be  used  in  the  analysis  of  fluid  transmission 
lines  by  considering  pressure  and  volumetric  flow  rate 
in  fluid  systems  to  be  analogous  to  potential  and  current 
respectively  in  electrical  systems.  With  this  analogy 
Ohms  law  becomes 


or 


where 


AP 

Q 

z 

Y 


AP  =  ZQ 

Q  =  YAP 

=  change  in  pressure 
=  volumetric  flow  rate 
=  complex  series  impedance 
=  complex  shunt  admittance 


(1) 


(2) 


=  1/Z 


The  complex  impedance  is  composed  of  the  resistance 
and  inertance,  while  the  admittance  is  composed  of  the  con¬ 
ductance  and  capacitance.  Using  these  components,  the 


v-  V 


above  equations  may  be  expressed  in  terms  of  the  distance 
along  the  line  as 


dx 


ZQ 


(R  +  jtoDQ 


(3) 


and 


(G  +  jwOP 


(4) 


where 

R  =  resistance 
L  =  inertance 
G  =  conductance 
C  =  capacitance 
all  taken  per  unit  length. 

Based  on  these  equations,  the  analogous  electrical 
circuit  to  a  fluid  transmission  line  is  shown  in  Figure  1. 


Blocked  and  Laminar  Flow 

The  equations  developed  by  Nichols  (Ref  1)  provide 
a  means  of  calculating  values  of  Z  and  Y  for  circular  lines 
with  either  laminar  mean  flow  or  blocked  line  conditions. 
These  equations,  however,  involve  the  use  of  Bessel  func¬ 
tions,  and  the  calculations  can  become  complicated.  Krish- 
naiyer  and  Lechner  (Ref  2)  have  presented  algebraically 
simplified  equations  which  provide  accurate  results  over 
the  frequency  range  0.1  <  w  <  <». 

frequency  ui  ,  is  defined  as 


The  characteristic 


The  transmission  line  may  now  be  described  by  the 
propagation  operator  T  and  the  characteristic  impedance  Z 


r  =  (ZY)^  =  ct  +  jB  (14) 

Z 

Zc  =  <15) 

where 

a  =  attenuation  per  unit  length 
S  =  phase  shift  per  unit  length 


For  the  general  case  of  a  line  with  terminal 
impedance  Z  ,  as  shown  in  Figure  2,  Kirshner  and  Katz 
(Ref  7)  express  the  input  impedance  as 


Z  . 


in 


Tz, 


cosh  (r£) 
sinh  (r£) 


+  Z^  sinh(r£) 

+  z^  cosh(re)j 


(16) 


The  pressure  transfer  function  relating  the  pres¬ 
sure  at  a  distance  x  along  the  line  to  the  input  pressure 
is 


„  sinhr  f,  (1-x/ £)  +  (Z^ /Z  )  coshF  ?,  (1-x/ 7 ) 
r*  _  _ L  C _ 

P.  sinhr?,  +  (Z^ /Z  )  coshm 

in  L  c 


(17) 


and  the  pressure  gain  is  defined  as  the  magnitude  of  the 


transfer  function 


or,  in  decibels. 


^dB  "  ^ 


(19) 


The  phase  shift  is  defined  as  the  angle  formed  by 
the  real  part  and  the  complex  part  of  the  transfer  function 
in  the  complex  plane 


B 


tan 


-1 


I  (P/P  ) 
m  in 


Re (P/P.  ) 
in 


(20) 


Turbulent  Flow 

While  Equations  (6)  and  (7)  have  proven  accurate 
for  use  in  blocked  line  and  laminar  flow  cases,  it  has  been 
shown  that  their  use  in  the  turbulent  flow  case  results 
in  attenuation  values  lower  than  those  observed  experi¬ 
mentally  (Refs  3-5) . 

The  velocity  profile  for  fully  developed  laminar 
steady  state  flow  is  parabolic  with  the  boundary  layer 
thickness  equal  to  the  tube  radius.  When  a  periodic 
signal  is  applied  to  the  steady  state  flow  a  dynamic  boun¬ 
dary  layer  is  created.  As  the  frequency  of  the  signal  is 
increased  this  dynamic  boundary  layer  becomes  thinner  and 
the  velocity  gradients  near  the  wall  increase.  As  a  result 
the  dynamic  layer  is  thinner  than  the  steady  state  boundary 
layer  at  frequencies  of  practical  interest  and  dictates 
the  line  response  (Ref  13) . 


In  turbulent  flow  the  steady  state  velocity  pro¬ 
file  exhibits  a  dependency  on  the  Reynolds  number.  The 
boyndary  layer  thickness  decreases  with  increasing  Reynolds 
number  and  at  the  limit  of  an  infinite  Reynolds  number  a 
uniform  profile,  or  slug  flow,  is  obtained.  With  this 
variation  in  steady  state  boundary  layer  thickness,  the 
dynamic  boundary  layer  created  by  the  application  of  a 
periodic  signal  will  no  longer  necessarily  be  thinner  than 
the  steady  state  boundary  layer.  With  the  response  deter¬ 
mined  by  the  thinner  of  the  two  layers,  it  is  clear  that 
variations  in  Reynolds  number  and  signal  frequency  would 
affect  line  response. 

In  order  to  classify  the  response  based  on  Reynolds 
number  and  frequency.  Brown  et  al.  (Ref  3)  presented  a 
plot  of  nondimensional  attenuation  versus  nondimensional 
frequency.  An  example  of  this  type  of  plot  is  shown  in 
Figure  3.  The  two  nondimensional  parameters  are  defined  as 
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nondimensional  attenuation  =  -  (21) 

V 

2 

nondimensional  frequency  =  -  =  8 —  (22) 

u 

This  plot  has  been  divided  into  three  frequency  regions 
by  Mooie  (Ref  4)  as  suggested  by  Brown  et  al.  (Ref  3) . 

In  the  low  frequency  region,  the  dynamic  boundary  layer  is 
thicker  than  the  steady  state  boundary  layer.  In  the 
mid  frequency  region,  the  dynamic  layer  and  the  steady 


layer  are  of  comparable  thickness.  In  the  high  frequency 
region  the  dynsimic  layer  has  become  much  thinner  than  the 
steady  state  layer  and  as  a  result  the  behavior  in  this 
region  is  the  same  as  in  laminar  flow.  The  frequency  at 
which  the  steady  state  layer  becomes  the  thicker  of  the 
two  and  the  response  transitions  to  laminar  behavior  is 
termed  the  break  frequency, 

O 

Several  turbulent  models  have  been  suggested  to 
predict  dynamic  response.  Moore  applied  a  constant  LRC 
model  to  the  low  and  mid  frequency  regions,  the  components 
of  which  were  selected  based  on  analysis  of  the  low  fre¬ 
quency  laminar  flow  case  (Ref  4) .  The  adiabatic  inertance 
was  defined  as 

=  p/A  (23) 
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A  steady  state  (DC)  resistance  was  defined  as 


DC 
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(24) 


and  using  the  Blasius  resistance  formula  for  turbulent  flow 
in  smooth  lines 
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This  was  rewritten  as 


The  capacitance  was  defined  as 


C  =  ^  (27) 

n  a 

where  n  is  the  polytropic  exponent  and  the  adiabatic 
capacitance.  C  is  defined  as 

3i 

C  =  A/kP  (28) 


The  value  of  n  is  seen  to  vary  from  k  at  high 
frequencies  to  1.0  at  low  frequencies. 

Wilkins  presented  modifications  to  the  constant 
LRC  model  in  an  attempt  to  improve  the  attenuation  predic¬ 
tions  of  his  model.  The  attenuation  curve  in  the  mid 
frequency  region  of  Figure  3  is  seen  to  be  dependent  only 
on  the  Reynolds  number.  The  modifications  to  the  resist¬ 
ance,  inertance  and  capacitance  proposed  by  Wilkins  were 
based  on  this  variation  in  attenuation  with  Reynolds  number 
A  linearized  (AC)  resistance  suggested  in  Ref  8 


AC 
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was  used  which  using  the  R^^^  in  Equation  (2  6)  reduces  to 


AC 


=  1.75  R 


DC 


(30) 


Moore  and  Franke  have  proposed  a  turbulent 
inertance  to  reflect  the  effect  of  the  velocity  profile  on 


the  inertance  (Ref  9) .  This  turbulent  inertance  is  defined 
as 


(31) 


(32) 


is  now  seen  to  be  a  function  of  the  Reynolds  number 
through  the  friction  factor  f. 

The  turbulent  capacitance  was  formed  by  combining 
the  imaginary  part  of  the  conductance  with  the  isothermal 
capacitance.  The  conductance  accounts  for  the  effects  of 
radial  heat  transfer  in  a  gas.  In  laminar  flow  the 
increased  attenuation  due  to  this  effect  is  well  under¬ 
stood  and  has  been  included  in  the  equations  by  Krishnaiyer 
and  Lechner  (Equations  (6)  and  (7)). 

The  conductance  for  turbulent  flow  as  presented 
by  Brown  et  al.  (Ref  3)  as 


{k-DC  L 


(33) 


If  the  quantity  is  assumed  to  be  much  less  than  L^, 

Wilkins  (Ref  6)  showed  that  after  some  algebraic  manipula¬ 
tion,  this  reduces  to 


G  =  ik-1) 
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From  this  expression,  the  turbulent  capacitance 
is  written  as 
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which  is  again  seen  to  be  dependent  on  the  Reynolds  number 
through  the  inertance  constant. 

The  second  turbulent  model  presented  by  Wilkins 
(Ref  6)  was  the  constant  LRCG  model.  This  model  uses  the 
same  resistance,  inertance  and  capacitance  components  as 
the  modified  constant  LRC  but  also  includes  the  real  part 
of  the  conductance  from  Equation  (34) .  Each  component  of 
these  models  is  seen  to  be  independent  of  the  signal  fre¬ 
quency  and  a  function  of  only  the  Reynolds  number. 

AC 

If  the  quantity  cannot  be  assumed  to  be  much 

less  than  LT,  the  complete  equation  for  G  must  be  utilized 
(Equation  (33) )  .  Expansion  of  this  expression  gives 


It  should  be  noted  that  use  of  the  full  expression 


results  in  a  component  that  is  dependent  on  frequency  as 
well  as  the  Reynolds  number. 

If,  as  before,  a  portion  of  the  imaginary  part  of 
G  is  combined  with  to  form  the  turbulent  capacitance, 
the  following  result  is  obtained. 
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The  turbulent  capacitance  thus  formed  is  also  seen  to  be 
frequency  dependent  in  contrast  to  the  component  used  by 
Wilkins  in  the  Constant  LRCG  model. 

A  new  model  can  be  formed  by  using  the  frequency 
dependent  expression  for  G  and  in  place  of  the  respec¬ 

tive  frequency  independent  terms  in  the  Constant  LRCG 
model.  The  resulting  expressions  will  be  referred  to  as 
the  Modified  LRCG  model. 

It  must  be  noted  that  although  the  expression  for 
the  conductance  used  in  both  LRCG  models  is  based  upon 
that  presented  by  Brown  et  al .  (Ref  3),  there  is  an 


important  difference  between  the  components  as  defined  by 
Brown  and  those  used  in  this  study.  As  an  example.  Brown 
defines  the  effective  inertia,  I  (equivalent  to  the 
adiabatic  inductance)  as 


I  =  P 


(38) 


as  compared  to  the  definition  of  L  in  Equation  (23) .  The 


two  components  are  seen  to  differ  by  the  area  factor  such 
that 


I  =  L  A 
a 


(39) 


A  similar  multiplicative  factor  must  be  included  in  the 
resistance  components  used  in  this  study  to  obtain  Brown's 
results . 

Similarly,  Brown's  admittance  components  are 
related  to  those  used  in  this  study  by 


Brown 


(40) 


Though  the  admittance  and  impedance  terms  each 
differ  by  a  factor  of  the  area,  these  factors  cancel  out 
when  the  two  components  are  combined  to  form  the  propaga¬ 
tion  operator  r  producing  the  same  result  regardless  of 
the  method  that  is  used.  Care  must  be  taken,  however,  to 
ensure  that  the  two  types  of  components  are  not  combined 
without  regard  to  the  multiplicative  factor  by  which  they 
differ . 


Terminal  Impedance 


Several  methods  have  been  proposed  to  accurately 
model  the  terminal  impedance  of  a  dynamic  system.  Using 
a  linearized  DC  expression  such  as  that  proposed  by  Franke 
Malanowski  and  Martin  (Ref  10) 


(41) 


shows  that  in  the  limiting  cases,  for  an  open  line, 

APt  =  0  and  Z  =0;  while  for  the  blocked  line  case, 

Li  Li 

AQt  =  0  and  Z^.  =  “>. 

Li  Li 

In  this  study  two  primary  methods  have  been  used 
to  model  terminal  impedance.  In  the  open  line  case,  the 
radiation  impedance  as  suggested  by  acoustic  theory  has 
been  employed. 

The  air  at  the  end  of  the  transmission  line  is 


described  as  a  massless  piston  which  transmits  some  of 
the  energy  into  free  space  while  reflecting  some  energy 
back  into  the  line.  This  finite  impedance  is  seen  to 
differ  from  the  zero  impedance  resulting  from  the  linear¬ 
ized  impedance  model  mentioned  above.  Morse  and  Ingard 


where  0  and  y  ^re  radiation  functions.  For  the  case  of 
o  ^o 

a  flanged  exit  (a  flush  opening  in  an  infinite  wall)  these 
functions  may  be  approximated  as 

9=0  (43) 

o 

X=0.4D(jj/a  (44) 

o 


Using  these  approximate  values  results  in  an  expression 
for  the  terminal  impedance  of 

=  -0.4j(jopD/A  (45) 

The  second  method  used  to  model  the  terminal 
impedance  involves  the  use  of  reflection  coefficients. 

The  orifice  of  the  transmission  line  was  modeled  as  illus¬ 
trated  at  Figure  4.  Kensler  and  Frey  (Ref  12)  give  the 
terminal  impedance  for  this  situation  as 
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where  W  is  the  reflection  coefficient  for  the  sudden 
change  in  line  area  between  1  and  2.  The  reflection 
coefficient  for  the  situation  where  line  2  is  terminated 
by  its  characteristic  impedance  is  given  as 

W  =  (Aj^-A^)  /  (A^+A2)  (47) 
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Schlichting  (Ref  13)  gives  the  presssure  drop  along  a 
line  as 

AP  =  -  f(Ax)  p  (u^/2D)  (49) 

and  by  substituting  the  applicable  expression  for  the 
friction  factor,  the  following  results  may  be  obtained. 

laminar  flow  AP  =  -0.106  (Ax)  pM/D^p  (50) 

turbulent  flow  AP  =  -(7.18x10  ^)  Ax  u ' * ^^p 

(51) 

From  these  expressions  it  can  be  seen  that,  for  a 
given  fluid  at  ambient  conditions  and  a  given  line  con¬ 
figuration,  the  mass  flow  rate  is  dependent  only  on  the 
pressure  drop  along  the  line.  In  test  situations  where 
an  upstream  bypass  device  is  used  to  simulate  a  branch 
line  in  a  more  complex  system,  the  pressure  drop  along  the 
test  line  may  be  used  to  establish  a  required  mass  flow 
rate. 

Error  Analysis 

In  order  to  properly  assess  the  correlation  of 
theoretical  predictions  with  experimental  results,  a  deter¬ 
mination  of  the  error  due  to  the  readability  of  measurement 
instruments  is  necessary.  In  this  study  the  method  of 
absolute  error  was  used  to  calculate  this  error  (Ref  15) . 

For  some  quantity  F,  defined  by 

F  —  F(ej^,e2,.-.,e^)  (52) 
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where  e^,  •••»  are  independent  measured  quantities, 

the  absolute  error  in  F  due  to  instrviment  readability  is 


AF 


(53) 


where  Ae^  is  the  readability  of  the  instriiment  used  to 

measure  e . . 

1 

The  absolute  error  may  be  used  to  determine  the 
accuracy  of  experimental  results.  If  the  error  analysis 
shows  an  error  of  5  percent  in  a  calculation,  for  example, 
the  effort  spent  on  resolving  a  difference  of  2  percent 
between  experimental  and  theoretical  results  might  be 
better  spent  improving  instrumentation  or  experimental 
techniques . 


Computer  Simulation 

The  computer  program  developed  by  Malanowski 
(Ref  14)  and  modified  by  Briski  (Ref  5)  and  Wilkins 
(Ref  6)  was  used  to  provide  theoretical  simulations  for 
comparison  with  experimental  results. 

The  test  line  is  divided  into  several  short  seg¬ 
ments  over  which  the  properties  are  assumed  to  remain 
constant.  Using  Schlichting ' s  equation  for  the  pressure 
drop  along  a  line  (Equation  (49) )  and  the  respective 
expressions  for  the  friction  factor  for  laminar  and  turbu¬ 
lent  flow,  Wilkins  developed  expressions  to  determine  the 
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pressure  drop  along  a  line  segment  of  known  length  given 
the  mass  flow  rate  and  the  inlet  pressure.  The  final 
expressions  are,  for  laminar  and  turbulent  flow  respec¬ 
tively. 


-  -  ReT  T 

AP  =  P.  -  P.  +64xu  - 

in  in 


P  -  Re^'^^T  T~l^ 

AP  =  P.  -  P.  ^+0.3164xu  - 5-^ 

in  in  jj3 


Using  the  appropriate  expression  the  program  calcu 
lates  the  pressure  drop  for  each  segment  in  succession 
and  thereby  the  inlet  pressure  for  each  of  the  line  seg¬ 
ments.  The  program  then  uses  the  segment  pressures  and 
other  properties  to  calculate  the  propagation  operator  and 
from  this,  the  theoretical  gain  and  phase  shift. 

The  program  is  interactive  and  offers  the  choice 
between  the  laminar  model  (Equations  (6)  and  (7))  and 
either  the  Modified  Constant  LRC  or  Modified  LRCG  models. 
For  the  turbulent  models,  the  attenuation  is  calculated 
by  both  the  selected  turbulent  model  and  the  laminar  model 
The  program  compares  the  two  values  and  selects  the  larger 
for  the  calculation  of  gain  and  phase  shift.  In  this 
manner  the  program  is  able  to  identify  the  break  frequency 
between  the  mid  frequency  and  high  frequency  ranges  and 
then  use  the  laminar  model  for  calculations  beyond  the 


break  frequency. 


The  program  also  contains  an  interactive  module 
to  allow  the  selection  of  the  appropriate  reflection 
coefficient  expression.  The  radiation  impedance  is  used 
for  open  line  cases. 

The  expressions  for  the  ratio  of  turbulent/ 
adiabatic  capacitances  and  the  inertance  constant  were 
added  to  the  program  to  reduce  the  amount  of  input  data 
required. 

The  documentation  and  layout  of  the  program  were 
improved  to  make  it  more  easily  understood  and  to  provide 
a  more  useful  analytical  tool  for  future  studies  in  this 
area.  A  complete  listing  of  the  program  is  included  at 
Appendix  B. 


Ill .  Experiment 


Apparatus 

The  test  sections  used  in  this  study  consisted  of 
24-inch  lengths  of  stainless  steel  line  which  were  per¬ 
manently  mounted  in  plexiglass  blocks  at  each  end.  Three 
inside  line  diameters  were  used:  0.041,  0.119  and  0.195  in. 
The  mounting  blocks  were  machined  to  provide  fittings  for 
air  input  line,  static  pressure  transducer,  and  dynamic 
pressure  transducer.  The  bypass  device  and  several  features 
of  the  upstream  mounting  block  as  listed  below  are  shown  in 
Figure  5. 

A  -  Bypass  device 
B  -  Static  pressure  transducer 
C  -  Dynamic  pressure  transducer 
D  -  Test  line 

A  typical  test  line  is  shown  in  Figure  6  and  the  actual  line 
dimensions  are  listed  in  Table  I.  The  downstream  mounting 
block  was  also  machined  to  allow  the  fitting  of  different 
termination  devices.  There  were  four  termination  devices 
for  each  test  line;  a  blocked  line  fitting,  an  open  line 
fitting  and  two  orifice  fittings  with  diameters  of  0.0135  in 
and  0.016  in.  (Figure  7).  Dimensions  of  the  termination 
devices  are  also  shown  in  Table  I. 

A  plexiglass  calibration  block  shown  in  Figure  8 


was  used  to  calibrate  the  dynamic  pressure  transducer/ 


ORIFICE  DEVICE  DIMENSIONS 


charge  amplifier  combinations.  The  calibration  block 
positioned  the  transducers  at  the  same  point  in  the  line 
to  allow  comparison  of  the  measured  signals.  The  block 
was  machined  similar  to  the  test  line  mounting  blocks  for 
an  air  supply  line  connection. 

The  experimental  instrumentation  used  in  the  study 
is  shown  in  Figures  9  and  10.  The  equipment  is  also  illus¬ 
trated  in  schematic  form  in  Figure  11. 

The  sending  wave  analyzer  provided  a  sinusoidal 
input  signal  of  variable  frequency,  which  was  measured  by 
a  frequency  counter,  to  the  amplif ier/power  supply  and  the 
pneumatic  driver.  This  input  signal  was  suitably  condi¬ 
tioned  and  passed  to  the  pneumatic  driver  itself  which 
imposed  a  sinusoidal  pressure  signal  on  the  mean  air  flow 
passing  through  the  driver.  The  pressure  signal  was  mea¬ 
sured  at  upstream  and  downstream  locations  by  piezoelectric 
dynamic  pressure  transducers.  The  charge  output  of  each 
dynamic  transducer  was  fed  to  a  charge  amplifier  where  it 
was  converted  to  a  voltage  and  input  to  a  wave  analyzer 
where  it  was  displayed  as  a  RMS  voltage.  The  charge  ampli¬ 
fier  output  was  also  input  to  a  digital  recording  oscillo¬ 
scope  where  the  waveforms  could  be  compared  and  recorded 
if  desired.  A  static  pressure  transducer  mounted  in  the 
upstream  mounting  block  provided  a  measure  of  the  static 
line  pressure  at  that  point.  The  pneumatic  driver,  test 
line  and  charge  amplifiers  are  shown  in  Figure  12. 
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Fig.  10.  Nicolet  Digital  Oscilloscope 


For  those  runs  involving  bypass  flow  upstream  of 
the  test  line,  a  complete  bypass  device  (Figure  5)  was 
used.  This  device  consisted  of  a  tee  fitting  located 
between  the  pneumatic  driver  and  the  upstream  mounting 
block.  The  bypass  outlet  of  the  tee  was  fitted  with  an 
inside  threaded  cap  with  a  floating  disc.  Adjustment  of 
the  threaded  cap  allowed  varying  amounts  of  the  flow  to 
escape  providing  increased  flow  through  the  pneumatic 
driver  for  a  given  flow  through  the  test  line.  This  device 
was  also  used  for  the  blocked  line  cases  where  the  entire 
flow  was  diverted  from  the  test  line. 

A  calibrated  rotometer  was  used  to  measure  the  mass 
flow  rate  through  the  test  line.  As  the  presence  of  the 
rotometer  downstreeim  of  the  test  line  altered  the  line's 
response,  flow  measurements  were  made  and  the  rotometer 
removed  before  data  was  taken.  When  the  bypass  device  was 
being  used,  the  rotometer  again  could  not  be  used  as  its 
presence  changed  the  respective  flows  through  the  test  line 
and  the  bypass  device.  To  determine  the  mass  flow  through 
the  test  line  with  the  bypass  device  open,  it  was  assumed 
that  the  difference  between  a  given  upstream  static  plea¬ 
sure  and  atmospheric  pressure  resulted  in  a  given  flow 
through  the  line. 

With  the  bypass  device  closed,  the  air  supply 
was  adjusted  to  give  the  desired  static  pressure  reading. 
The  rotometer  was  then  used  to  measure  the  flow  and 
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removed.  The  bypass  device  was  then  opened  and  the  air 
supply  increased  until  the  previous  static  pressure  read¬ 
ing  was  re-established. 

A  complete  listing  of  the  instrxamentation  and 
ancillary  equipment  used  in  this  study  is  provided  in 
Appendix  C. 

Procedures 

Dynamic  Transducer  Calibration.  Throughout  the 
study  each  dynamic  pressure  transducer  was  paired  with  a 
specific  charge  amplifier  and  the  combination  treated  as  a 
single  unit.  To  calibrate  the  transducer/amplif ier  units 
the  transducers  were  mounted  in  the  calibration  block  which 
was  then  connected  to  the  air  supply  line  in  place  of  the 
test  line.  The  remainder  of  the  experimental  equipment 
was  set  up  as  for  a  normal  test  run.  At  a  selected  fre¬ 
quency,  the  signals  from  the  transducers  were  displayed 
on  the  oscilloscope  and  positioned  such  that  the  signals 
overlayed  each  other.  With  the  magnitude  switches  for 
each  channel  set  to  the  same  setting,  the  sensitivity  dial 
on  the  downstream  charge  amplifier  was  adjusted  until  the 
signal  traces  on  the  oscilloscope  exactly  matched.  The 
amplitude  of  the  signals  was  verified  by  setting  the  down¬ 
stream  wave  analyzer  frequency  to  obtain  maximum  voltage. 


It  was  found  that  using  the  procedure  of  matching  the 
signals  visually  on  the  oscilloscope  produced  voltages  on 
the  wave  analyzers  which  were  equal  to  within  the  accuracy 
of  the  analyzer  meters. 

This  sensitivity  setting  was  then  verified  over  a 
range  of  frequencies  and  air  supply  levels.  In  all  cases 
the  signals  displayed  negligible  differences  over  the 
range  of  frequency  and  air  supply  levels  used  in  this 
study.  After  calibration  the  transducers  were  transferred 
to  the  applicable  test  line  and  the  line  connected  in  place 
of  the  calibration  block. 

This  calibration  procedure  was  repeated  each  time 
the  transducers  were  placed  in  a  new  test  line.  Initial 
tests  using  the  calibration  block  showed  that  the  same 
results  were  obtained  regardless  of  the  position  in  which 
the  transducers  were  placed.  Calibration  tests  were 
repeated  with  the  transducer  positions  reversed  without 
measurable  differences. 

Experimental  Procedure.  The  experimental  data 
collected  during  this  study  consisted  of  the  sending 
(upstream)  pressure,  the  receiving  (downstream)  pressure 
and  the  phase  shift  between  the  two  pressure  signals.  The 
data  was  measured  as  a  function  of  the  frequency  over  the 
frequency  range  20  to  1000  Hz. 


Given  a  selected  Reynolds  number,  the  required 
static  pressure  setting  was  determined  from  Equations  (50) 
or  (51) .  Following  the  warmup  period  of  the  pneumatic 
driver  amplifier /power  supply  and  the  zeroing  of  the  wave 
analyzers,  the  air  pressure  was  adjusted  to  obtain  the 
correct  reading  on  the  voltmeter  connected  to  the  static 
pressure  transducer.  The  voltage  reading  was  converted 
to  psig  using  the  calibration  equation  determined  for  the 
transducer.  The  ambient-atmospheric  pressure  and  tempera¬ 
ture  were  recorded  with  the  atmospheric  pressure  being 
corrected  for  temperature  and  gravity  effects. 

The  frequency  was  selected  on  the  sending  wave 
analyzer  using  the  frequency  counter  as  the  measuring 
device.  The  frequency  on  the  receiving  wave  analyzer  was 
adjusted  until  a  maximum  voltage  was  obtained.  This  pro¬ 
cedure  was  used  rather  than  simply  matching  the  frequen¬ 
cies  using  the  frequency  counters  because  it  was  dis¬ 
covered  early  in  the  study  that  there  were  sometimes  differ 
ences  of  up  to  7  Hz  between  the  sending  frequency  and  the 
receiving  frequency  at  which  the  maximum  voltage  was  mea¬ 
sured.  In  many  instances  if  the  frequencies  were  simply 
matched,  the  receiving  pressure  signal  would  decrease  from 
a  maximum  of  anywhere  from  50  to  300  mV  to  almost  zero  mV. 
This  problem  is  attributed  to  the  wave  analyzers  used  in 
this  study.  The  problem  should  be  eliminated  by  the 
scheduled  acquisition  of  new  wave  analyzers. 


When  the  maximum  voltage  was  obtained,  both  charge 
amplifiers  were  grounded  using  the  ground  button  on  the 
front  panel  to  ensure  accurate  readings  and  the  two  pressure 
signal  voltages  recorded. 

The  storage  feature  of  the  Nicolet  oscilloscope 
was  utilized  to  record  the  traces  of  the  pressure  signals. 
The  time  difference  between  identifiable  common  points  on 
the  two  signals  was  determined  using  the  movable  cursor 
and  the  digital  display  of  the  oscilloscope.  Any  traces 
which  showed  distinctive  or  unusual  features  were  stored 
on  the  oscilloscope's  magnetic  disc  for  later  review  and 
analysis.  This  procedure  was  repeated  over  the  usable 
frequency  range  of  the  pnexmatic  driver  of  20  to  1000  Hz 
at  varying  intervals. 

The  flow  rate  through  the  test  line  was  measured 
using  the  rotometer  before  and  after  the  test  run.  The 
rotometer  measurement  was  converted  to  a  mass  flow  rate 
using  calibration  curves  previously  prepared  for  the 
rotometer.  The  calibration  of  the  most  commonly  used 
rotometer  tube  was  verified  during  this  study. 

The  air  supply  was  monitored  using  the  static 
pressure  transducer  output.  Fluctuations  in  the  pressure 
signal  voltage  readings  could  be  related  to  sharp  fluctua¬ 
tions  in  the  supply  line  pressure.  Data  measurements  were 
not  taken  during  periods  of  line  pressure  fluctuation. 


The  precaution  of  grounding  the  charge  amplifiers 
before  taking  the  readings  was  implemented  after  it  was 
noted  that  the  output  of  the  amplifiers  would  decay  sig¬ 
nificantly  over  a  period  of  time.  Although  the  length  of 
time  required  for  this  effect  to  occur  was  longer  than 
that  usually  taken  to  obtain  the  readings,  this  precaution 
was  taken  to  ensure  accurate  readings  in  the  event  that 
there  was  some  delay. 

In  the  cases  where  the  bypass  device  was  used,  the 
required  static  pressure  reading  was  obtained  with  the 
bypass  device  closed  by  adjusting  the  air  supply  and  the 
flow  rate  was  measured.  The  bypass  device  was  then  opened 
the  optimum  amount  to  provide  maximum  pressure  signal 
voltages  and  clear  oscilloscope  traces  and  the  air  supply 
increased  to  re-establish  the  static  pressure  reading. 

Test  Designation 

A  five-digit  numeric  code  was  assigned  to  each 
test  run.  The  code  was  the  same  as  that  used  by  Wilkins 
(Ref  6)  to  allow  direct  comparison  of  experimental  results 
and  to  maintain  consistency  within  the  data  base. 

The  first  digit  identifies  the  line  used;  1  for 
the  0.041,  2  for  the  0.119  and  3  for  the  0.195  in.  diameter 
line.  The  second  digit  identifies  the  line  termination 
condition;  0  for  blocked,  1  and  2  respectively  for  the 


0.0135  and  0.016  in.  diameter  orifices,  and  3  for  open. 

The  last  three  digits  indicate  the  nominal  Reynolds  number 
in  hundreds. 

An  alpha  character  suffix  was  used  in  some 
instances  to  identify  significant  test  characteristics. 

The  characters  and  their  meaning  are  listed  below, 
b  -  bypass  device  open 
n  -  modified  orifice  device 
r  -  transducer  positions  reversed 


IV.  Results  and  Discussions 


Test  Conditions 

A  sxairanary  of  line  configurations  and  ambient  con¬ 
ditions  for  which  experimental  tests  were  conducted  is 
presented  in  Table  II .  The  shaded  area  of  Figure  13  shows 
the  portion  of  the  nondimensional  frequency  vs.  attenua¬ 
tion  region  covered  by  these  results  (3.52  <  (^  <  9201.8). 

Analysis  of  Instrument  Error 

Table  III  presents  the  absolute  error  in  the  experi¬ 
mental  gain  for  typical  pressure  measurements.  This  error 
is  based  solely  on  the  readability  of  the  wave  analyzers 
and  was  determined  by  analysis  of  Equation  (19)  using  the 
method  of  Equation  (53). 

A  similar  analysis  was  conducted  for  the  phase 
angle  calculation  based  on  the  readability  of  the  Nicolet 
oscilloscope  and  the  frequency  counters.  The  experimental 
phase  shift  is  calculated  from  the  equation 

B  =  -0.36 (freq) (phd)  (5G) 

where  freq  =  signal  frequency  in  Hz 
phd  =  phase  delay  in  msec. 

The  resultant  absolute  error  for  typical  phase  delay  mea¬ 
surements  is  presented  in  Table  IV. 
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TABLE  III 


ABSOLUTE  ERROR  IN  EXPERIMENTAL  GAIN 
DUE  TO  INSTRUMENT  READABILITY 


Error 


2 

22.0 

11.0 

.34 

3.1 

9 

27.0 

-8.1 

.29 

3.5 

0 

410 

5.02 

.40 

8.0 

0 

340 

-6.64 

.37 

5.6 

5 

7.2 

-2.4 

.21 

8.8 

6 

81 

4.9 

.30 

6.0 

0 

200 

-3.52 

.36 

10.3 

870 

-5.82 

.36 

6.1 

5 

1450 

-1.122 

.56 

50.2 

Freq 

(Hz) 


TABLE  IV 

ABSOLUTE  ERROR  IN  EXPERIMENTAL  PHASE  SHIFT 
DUE  TO  INSTRUMENT  READABILITY 


(msec) 


1.2 

2.5 

2.0 

0.58 

0.78 


Phase  Shift 
(deg) 


-25.9 

-90.0 

-331.2 

-510.3 

-601.5 


Error 


(deg) 


21.3 

6.8 

1.6 

1.1 

1.2 


-■V 


' 


Though  the  Nicolet  oscilloscope  features  sample 
-9 

rates  up  to  50  x  10  secs.,  the  sample  rates  used  in  this 
study  were  determined  largely  by  the  signal  frequency. 

It  was  necessary  to  select  a  sample  rate  which  provided 
a  display  of  several  cycles  to  allow  a  significant  number 
of  phase  delay  measurements  to  be  made  at  each  frequency. 

The  major  source  of  error  in  the  phase  shift  cal¬ 
culation  is  the  determination  of  the  actual  peak  location 
on  the  upstream  signal  in  certain  situations.  In  these 
instances  the  signal  either  degenerated  into  a  multipeak 
signal  as  illustrated  in  Figure  14  or  displayed  random 
noise  which  obscured  the  location  of  the  peak.  The  range 
of  possible  peak  locations  in  these  situations  was  several 
sample  steps  wide  and  presents  significant  error  in  addi¬ 
tion  to  that  resulting  from  instrument  readability.  Almost 
all  cases  where  the  experimental  phase  shift  data  differs 
significantly  from  theoretical  predictions  can  be  traced 
to  situations  where  this  problem  was  evident. 

General  Gain  and  Phase 
Shift  Results 

All  cases  presented  in  Table  II  were  evaluated 
using  the  computer  simulation  to  compare  experimental 
results  with  theoretical  predictions.  Figures  15-32 
present  the  correlation  between  experimental  and  theoreti¬ 
cal  results  for  some  representative  cases  from  this  study. 
The  results  of  the  remaining  cases  will  be  used  to 
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illustrate  points  discussed  in  the  following  sections. 

In  most  cases  the  experimental  results  show  good  agreement 
with  the  results  of  Brisk!  (Ref  5)  and  Wilkins  (Ref  6) . 
Blocked  line  and  laminar  mean  flow  experimental  results 
show  close  correlation  with  laminar  theory  predictions 
(Figures  15-18  and  27-30) ,  while  in  turbulent  mean  flow 
cases  the  laminar  theory  predicts  significantly  lower 
attenuation  than  observed  experimentally  (Figures  19-26 
and  31-32) .  The  Modified  Constant  LRC  model  (MCLRC) 
developed  by  Wilkins  (Ref  6)  predicts  increased  attenua¬ 
tion  in  turbulent  mean  flow  cases  but  the  attenuation  is 
still  below  experimental  results.  The  use  of  the  Constant 
LRCG  model  (Ref  6)  results  in  a  further  increase  in  attenu¬ 
ation  over  that  predicted  by  the  MCLR  model,  but  in  most 
cases  the  theoretical  attenuation  still  results  in  peak 
gain  magnitudes  less  than  those  observed  experimentally. 

The  phase  shift  results  are  generally  in  close  agreement 
between  experimental  and  theoretical  models  with  some 
deviation  at  high  Reynolds  numbers.  The  deviations  are 
consistent  with  results  obtained  by  Wilkins  (Ref  6)  but 
are  of  unknown  cause  at  this  time. 

Analysis  of  Radiation  Impedance 

For  the  open  line  cases  in  this  study  the  radia¬ 
tion  impedance  was  used  in  the  computer  simulation  to 
determine  the  terminal  impedance.  A  comparison  of  the 
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frequency  range.  The  increase  was  so  small,  however,  that 
it  is  considered  insignificant  in  relation  to  the  differ¬ 
ences  between  experimental  and  theoretical  results.  The 
effect  of  the  radiation  impedance  is  seen  to  be  negligible 


Figures  40-48  show  the  comparison  of  experimental  gain 
between  the  bypass  and  no-bypass  conditions. 

The  experimental  results  match  very  closely  in 
most  cases  confirming  that  flow  variations  upstream  of  the 
test  section  have  no  effect  on  line  response. 

The  results  of  Case  21004  (Figure  45) ,  however, 
show  an  unexpected  positive  offset  of  the  gain  over  the 
entire  frequency  range  in  the  bypass  condition.  This  off¬ 
set  resulted  from  a  switch  of  the  transducer  locations  and 
indicated  that  there  was  a  difference  in  sensitivity 
between  the  two  transducer/charge  amplifier  pairs  despite 
the  same  settings  being  used  on  both  amplifiers.  This 
difference  was  confirmed  by  the  results  of  two  other  sets 
of  runs  where  the  transducer  positions  were  reversed. 

These  results  led  to  the  implementation  of  the  calibration 
procedure  discussed  in  Chapter  III.  Results  obtained  after 
calibration  of  the  transducers  show  no  such  offset.  The 
theoretical  gain  has  not  been  plotted  for  those  cases  con¬ 
ducted  before  transducer  calibration  as  the  difference  in 
sensitivity  caused  an  offset  of  both  the  bypass  and 
no-bypass  results  from  the  actual  values. 

Effect  of  Modified  Orifice  Devices 

Due  to  the  machining  methods  used  in  their  produc¬ 
tion,  the  original  orifice  devices  did  not  match  the  char¬ 
acteristics  of  the  orifice  model  shown  in  Figure  4.  The 
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orifice  devices  contained  a  significant  length  of  area 
reduction  before  the  actual  orifice  as  shown  by  the  dotted 
lines  of  Figure  7.  The  devices  were  modified  to  more 
closely  match  the  model  by  machining  out  the  inlet  passage 
to  a  constant  dicimeter.  The  dimensions  noted  on  Figure  7 
and  listed  in  Table  I  apply  to  the  modified  orifice  devices. 

The  gain  and  phase  shift  results  in  Figures  49-60 
show  close  correlation  with  results  of  similar  cases  pre¬ 
sented  by  Wilkins  with  the  exception  of  Case  12034  (Figures 
66-67) .  The  experimental  gain  data  from  this  case  shows 
significant  deviation  from  Wilkins'  results  in  similar 
cases  (12020  and  12050) .  The  results  show  higher  peak 
gains  over  the  frequency  range  20-500  Hz.  Particularly 
noticeable  is  the  peak  gain  of  6.89  dB  at  350  Hz  compared 
with  values  of  approximately  zero  in  Wilkins’  results. 

The  gain  results  over  the  remainder  of  the  frequency  range 
and  the  phase  shift  results  show  good  agreement  with 
Wilkins'  results.  The  deviation  between  the  experimental 
results  and  the  theoretical  predictions  is  generally  the 
same  as  displayed  in  Wilkins'  results.  The  reason  for  this 
anomaly  is  not  known  but  the  almost  exact  agreement  between 
Cases  12034  and  12034b  reduces  the  likelihood  of  experi¬ 
mental  error.  It  should  be  noted  that  there  is  a  signifi¬ 
cant  difference  in  the  Reynolds  numbers  between  this  case 
and  comparable  cases  of  Wilkins  (3425  vs.  2044  and  4992). 
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Modified  LRCG  Model 

Analysis  of  theoretical  results  obtained  using 
the  Constant  LRCG  model  developed  by  Wilkins  (Ref  6) 
revealed  that  this  model  did  not  display  the  characteris¬ 
tic  decrease  in  nondimensional  attenuation  at  low  non- 
dimensional  frequency  as  described  by  Brown  et  al.  (Ref  5)  . 
(All  further  references  to  attenuation  and  frequency  will 
refer  to  the  nondimensional  parameters  unless  otherwise 
specified. )  A  similar  analysis  of  the  MCLRC  model  showed 
that  the  MCLRC  model  did  exhibit  the  expected  decrease  in 
attenuation.  This  led  to  a  review  of  the  development  of 
the  Constant  LRCG  model  to  determine  the  cause  of  the 
deviation. 

The  assumption  that  RAC/uj  is  much  less  than  L^ 
was  used  by  Wilkins  to  simplify  the  model  and  remove  any 
dependency  on  the  frequency.  In  his  results  Wilkins 
(Ref  6,  p.  41)  discusses  the  limitations  imposed  on  the 
model  by  the  assiomption.  A  review  of  Wilkins'  cases 
reveals  that  for  a  wide  range  of  line  configurations  and 
Reynolds  numbers  this  assumption  is  not  valid  over  a  sig¬ 
nificant  portion  of  the  frequency  range  under  consideration 
Table  V  presents  a  review  of  cases  from  both  this  study 
and  Wilkins  and  illustrates  the  frequency  limitation  of 
this  assumption  for  various  line  configurations  and  flow 
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TABLE  V 


COMPARISON  OF  VS. 


Case 

Re 
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Ibf-s  4 
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lb  -s^  4 

,rad. 

(Hz) 

\  ^  ) xxu 
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in 

‘  s  ' 

11030* 

2968 

1434.8 

2.1388 

6283 

1000 

12018n 

1780 

906.96 

1.8190 

4986 

794 

23070* 

7326 

40  .542 

.10723 

3781 

602 

23090 

9094 

43.857 

.10892 

4027 

641 

23100* 

10450 

48.683 

.11197 

4348 

692 

23105 

10490 

48.812 

.11037 

4423 

704 

23165 

16513 

68.598 

.12297 

5578 

888 

31003b 

300 

.46812 

.086477 

54.1 

8.6 

31009 

927 

1.0746 

.11011 

97.6 

15.5 

31052 

5156 

3.9743 

.037176 

1069 

170 

32010n 

979 

1.1190 

.12665 

88.4 

14.1 

33020* 

1818 

1.8192 

.037723 

482 

76.8 

33070* 

7177 

5.5296 

.038131 

1450 

231 

33079 

7862 

5.3430 

.037187 

1437 

229 

*denotes  cases  from  Wilkins  (Ref  6). 


is  defined  as  the  frequency  at  which  RAC/w 
=  0.1 (L^)  and  above  which  the  assumption  RAC/w  much  less 
than  is  taken  to  be  valid. 


The  results  in  Table  V  show  that  the  assumption 
is  invalid  for  the  0.041  in.  diameter  line  over  most  of  the 
frequency  range  except  at  very  low  Reynolds  numbers.  In 
the  case  of  the  0.195  in.  diameter  line  the  assumption  is 
generally  valid  except  at  low  frequencies  (<200  Hz)  and 
high  Reynolds  numbers  (>7000) . 

Results  obtained  using  the  Modified  LRCG  model 
which  incorporates  the  full  expression  for  G  (Equation 
36))  display  the  expected  decrease  in  attenuation  with 
negligible  changes  in  attenuation  at  higher  frequencies. 
Figures  61-63  show  comparisons  of  the  attenuation  predicted 
by  the  Constant  LRCG  model  with  that  predicted  by  the  Modi¬ 
fied  LRCG  model.  The  difference  in  attenuation  at  low 
frequencies  is  immediately  obvious. 

Though  the  inclusion  of  the  frequency  dependent 
term  drastically  alters  the  nondimensional  behavior  of  the 
model,  the  effect  on  the  gain  and  phase  shift  predictions 
is  less  significant.  Figures  64-73  compare  the  effect 
of  the  Modified  LRCG  model  versus  the  Constant  LRCG  model 
for  typical  cases  from  this  study.  These  figures  indicate 
an  increase  in  the  gain  over  that  predicted  by  the  Con¬ 
stant  LRCG  model,  particularly  at  low  frequency.  The 
increase  in  the  gain  depends  on  the  line  configuration  and 
signal  frequency  and  generally  reflects  the  validity  of  the 
limiting  assumption  discussed  earlier.  Review  of  the  theo¬ 
retical  data  from  Case  11033  (Figure  64),  for  example. 


shows  an  increase  of  3.5  dB  at  1  Hz  and  of  0.9  dB  at  100  Hz. 


The  differences  decreases  with  increasing  frequency  and  in 

this  case  is  negligible  above  about  700  Hz.  This  increase 

in  the  gain  produces  better  agreement  with  experimental 

results  yet  a  significant  difference  still  remains  between 

theoretical  and  experimental  results. 

Examination  of  the  expression  for  the  modified 

(Equation  (37)  )  shows  that  as  oi  varies  from  0  to  infinity, 

k-1 

C„  varies  from  kC  to  (k  -  — )C  .  The  latter  value  is 

T  a  ^LT  ^ 

the  same  as  the  original  expression  for  (Equation  (35) ) . 
As  illustrated  in  Table  VI  the  use  of  Equation  (37)  results 
in  an  increase  in  of  from  approximately  38  percent  at 
1  Hz  to  less  than  1  percent  at  1000  Hz  over  the  value  of 
determined  by  Equation  (35) .  It  is  also  noted  that  the 
value  of  the  original  is  less  than  2  percent  greater 
than  C  for  the  cases  examined  in  this  study. 

Table  VII  shows  the  effect  of  the  modification  of 
on  the  attenuation  predicted  by  the  Modified  LRCG  model. 
To  obtain  this  data  the  Modified  LRCG  model  was  used  with 
values  of  calculated  from  Equations  (35)  and  (37)  . 

The  use  of  Equation  (37)  results  in  an  increase  in  the 
attenuation  over  that  predicted  using  Equation  (35) . 

Though  the  percentage  increase  is  large  at  low  frequencies, 
the  absolute  increase  is  on  the  order  of  2  dB  and  is  much 
less  significant  when  compared  with  the  25-55  dB  (300-600 
percent)  decrease  in  attenuation  resulting  from  the  use  of 
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TABLE  VI 

VARIATION  OF  Cm  DUE  TO  FREQUENCY 


(Eq  35) 

X  10^) 


Range  of  C  (Eq  37) 

in^  4 
(-^  X  10^)** 


11012 

1180 

.334 

.459  -  .337 

12018n 

1780 

.315 

.435  -  .318 

23070* 

7326 

5.285 

7.308  -  5.292 

23090 

9094 

5.151 

7.126  -  5.159 

23105 

10490 

5.038 

6.972  -  5.048 

23110* 

10450 

5.048 

6.986  -  5.058 

23165 

16513 

4.537 

6.286  -  4.551 

"denotes  case  from  Wilkins  (Ref  6) . 

**  (1  <  oj  <  1000)  . 

TABLE  VII 

PERCENT  INCREASE  IN  ATTENUATION  (MODIFIED  LRCG  MODEL) 
DUE  TO  MODIFIED  C„  (EQ  37)  OVER  C^  (EQ  35) 


Frequency 

(Hz) 


11033 


23090 


Case 

23105 


23165 


33079 


17.6 

17.5 

17.5 

17.6 

9.9 

7.8 

8.5 

10.4 

4.5 

3.0 

3.5 

4.9 

1.4 

0.84 

1.0 

1.6 

0.64 

0.38 

0.46 

0.73 

0.24 

0.14 

0.17 

0.27 

- 

0.04 

0.04 

0.07 

the  Modified  LRCG  model  vice  the  Constant  LRCG  model.  This 


indicates  that  the  largest  part  of  the  difference  in  attenu¬ 
ation  between  the  two  LRCG  models  is  due  to  the  use  of  the 
full  expression  for  Re  (G)  . 

Correlation  of  Modified  LRCG  Model 
with  Brown ' s  Modified  IRC  Model 

In  Reference  3,  Browr.  et  al.  present  the  results 
of  their  Modified  IRC  model  for  turbulent  flow  accounting 
for  heat  transfer.  The  laminar  reference  curve  (Figure  3) 
produced  using  the  simplified  equations  developed  by 
Krishnaiyer  and  Lechner  (Ref  3)  matches  the  laminar  refer¬ 
ences  shown  in  Figure  9  of  Reference  3  and  in  Wilkins 
(Ref  6)  except  for  the  difference  in  the  slope  of  the  refer¬ 
ence  compared  to  the  Modified  LRCG  curves  in  the  low  fre¬ 
quency  region.  Brown  shows  a  change  in  the  slope  of  the 
laminar  reference  at  =  2.5  which  does  not  appear  in  the 
laminar  reference  produced  using  the  equations  of 
Krishnaiyer  and  Lechner  (Equations  (6)  and  (7) ) . 

The  mid  frequency  range  attenuation  levels  predicted 
by  the  Modified  LRCG  model  agree  with  those  shown  on 
Brown's  figure  to  within  the  resolution  allowed  by  the 
figure . 

Brown  (Ref  3)  gives  an  expression  for  the  break 
frequency  between  the  mid  and  high  frequency  regions  as 

Q,  =  2.5  X  10“^  Re^'®  (57) 


This  expression  identifies  a  point  of  inflection  on 
Brown's  curves  rather  than  the  intersection  of  the  turbulent 
model  curves  and  the  laminar  reference  curve.  The  results 
of  this  study  give  the  expression  for  the  intersection 
point  for  the  Modified  LRCG  model  as 

=  9.68  X  10“'^  Re^*^®  (58) 

assiaming  a  power  form  solution  as  presented  by  Brown. 

This  break  frequency  agrees  closely  with  the  actual  inter¬ 
section  point  of  the  curves  in  Brown's  figure. 

Correlation  of  Results 
with  Nichols  (Ref  1) 

Figure  74  presents  a  comparison  of  the  attenuation 
predicted  by  the  modified  equations  of  Krishnaiyer  and 
Lechner  (Equations  (6)  and  (7))  with  that  given  by 
Nichols  (Ref  1,  Figure  6) .  The  large  deviation  shown  is 
attributed  to  the  use  of  an  approximation  to  the  Bessel 
functions  used  to  derive  Equations  (6)  and  (7) .  The  simpli¬ 
fied  equations  are  stated  to  be  limited  to  the  region 
Go/cj^  >  0.1.  This  is  confirmed  by  Figure  74. 

Figures  75  and  76  compare  the  attenuation  predic¬ 
tions  of  the  various  models  discussed  in  this  study.  The 
deviation  of  the  Constant  LRCG  at  low  frequencies  is  again 
seen  to  be  significant.  The  MCLRC  model  and  the  Modified 
LRCG  model  show  comparable  behavior  with  the  Modified  LRCG 
model  displaying  the  higher  attenuation  of  the  two  models. 


V.  Conclusions 


1.  Changes  in  conditions  upstream  of  the  test 
section  as  represented  by  a  bypass  device  with  variable 
flow  rates  did  not  affect  the  line  response  as  measured 
experimentally.  This  was  determined  to  be  true  in  both 
laminar  and  turbulent  mean  flow  conditions. 

2.  The  elimination  of  small  differences  between 
the  orifice  devices  and  the  orifice  model  used  to  determine 
the  reflection  coefficient  caused  no  change  in  the 
experimental  results. 

3.  The  absolute  error  in  the  experimental  calcu¬ 
lations  due  to  instrument  readability  was  determined  to 
be  approximately  0.5  dB  in  the  gain  calculations  and 
approximately  5  degrees  in  the  phase  shift  calculations. 
With  the  current  experimental  apparatus  these  values 
represent  realistic  limits  on  the  resolution  of  differences 
between  theoretical  predictions  and  experimental  results. 

4.  The  effect  of  the  radiation  impedance  on  the 
theoretical  predictions  in  open  line  cases  is  negligible. 
The  use  of  Z  =  0  rather  than  the  radiation  impedance 
resulted  in  an  insignificant  increase  in  the  gain. 

5.  The  Modified  LRCG  model,  which  is  formed  by 
adding  the  full  frequency  dependent  expression  for  the 


conductance  to  the  MCLRC  model,  was  seen  to  predict  lower 
attenuation  than  the  Constant  LRCG  model  at  low  frequency. 
The  resultant  increase  in  the  gain  varied  from  as  much  as 
6  dB  for  the  small  diameter  line  to  less  than  1  dB  for  the 
large  diameter  line.  The  change  was  negligible  at  high 
frequency. 
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Appendix  A:  Comparison  of  Experimental 
and  Theoretical  Results 

Case  Number  Explanation 
Case  ABxxx 
A  -  line  diameter 

1  =  0.041  in. 

2  =  0.119  in. 

3=0. 196  in. 

B  -  line  termination 
0  =  blocked 

1  =  0.0135  in.  diameter  orifice 

2  =  0.016  in.  diameter  orifice 

3  =  open 


XXX  =  Nominal  Reynolds  number  in  hundreds 
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Theoretical  and  Experimental  Phase  Shift — Case  10000  Laminar  Theory 
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FREQUENCY 


Fig.  26,  Theoretical  and  Experimental  Phase  Shift — Case  23165 
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Theoretical  and  Experimental  Gain — Case  31003b — Laminar  Theory 


Theoretical  and  Experimental  Phase  Shift — Case  31003b — Lciminar  Theory 


Effect  of  Radiation  Impedance  on  Theoretical  Gain — Case  23090 — Constant  LRCG 


Effect  of  Radiation  Impedance  on  Theoretical  Gain — Case  23165 — Constant  LRCG 


Effect  of  Radiation  Impedance  on  Theoretical  Gain — Case  33052 — Constant  LRCG 


Effect  of  Radiation  Impedance  on  Theoretical  Gain — Case  33060 — Constant  LRCG 


H/0  BYPRSS 


NOQ.  CONST. 


Effect  of  Upstream  Conditions  on  Experimental  Gain — Case  11014n 


H/0  BYPASS 


Effect  of  Upstream  Conditions  on  Experimental  Gain — Case  11020 


H/0  arPRSS 


Effect  of  Upstream  Conditions  on  Experimental  Gain — Case  11033 


HOD.  CONST 


TRflNbOUCERS 


HOO.  CONST. 


Effect  of  Upstream  Conditions  on  Experimental  Gain — Case  21007 


HOD.  CONST 


Effect  of  Upstream  Conditions  on  Experimental  Gain — Case  23107 


U/0  BTPnSS 


Effect  of  Upstream  Conditions  on  Experimental  Gain — Case  31009 


Theoretical  and  Experimental  Gain — Case  11030n 


FREQUENCir  (HERTZ)  «10‘ 

Theoretical  and  Experimental  Phase  Shift — Case  110 30n — Mod.  Const.  LRC 


Theoretical  and  Experimental  Gain — Case  120 18n — Laminar  Theory 


Theoretical  and  Experimental  Gain — Case  22020n — Laminar  Theory 
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Theoretical  and  Experimental  Phase  Shift — Case  22020n — Laminar  Theory 


Theoretical  and  Experimental  Phase  Shift — Case  31009n — Laminar  Theory 


Nondimensional  Behavior — Constant  LRCG  vs.  Modified  LRCG — Case  11033 


Nondimensional  Behavior — Constant  LRCG  vs.  Modified  LRCG — Case  33060 


Nondimen sional  B^ava^ — Constant  LRCG  vs 


Effect  of  Modified  LRCG  on  Theoretical  Gain — Case  11033 


Constant  LRCG 


Theoretical  Phase  Shift — Case  11033 


Effect  of  Modified  LRCG  on  Theoretical  Gain — Case  12034 
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Effect  of  Modified  LRCG  on  Theoretical  Phase  Shift — Case  23165 


Effect  of  Modified  LRCG  on  Theoretical  Phase  Shift — Case  33060 


Fig.  72.  Effect  of  Modified  LRCG  on  Theoretical  Gain — Case  33079 


Effect  of  Modified  LRCG  on  Theoretical  Phase  Shift — Case  33079 


AND  LECHNCR 


Fig.  74,  Comparison  of  Theoretical  Attenuation — Nichols  vs.  Krishnaiyer-Lechner 
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Fig.  76.  Comparison  of  Attenuation  for  Theoretical  Models — Case  33079 


the  Pneumatic  Transmission  Line 


PROGRAM  MSB 

C 

C  TRANSMISSION  LINE  STUDY/Z-IN  METHOD  WITH  BRANCH 

DIMENSION  DBT{300) ,RVT(100) ,DI (100) ,AD(100) ,AR(100) ,CV(100) , 
SOMTdOO)  ,CNA(100)  ,GMA(100)  ,AGM(100)  ,FN(100)  ,RN(100)  ,GN(100)  , 
SALNdOO)  ,CNdOO)  ,ANd00)  .BTNdOO)  .AMCdOO)  ,AZRNdOO)  ,BZBNd00)  , 
SAZINCIOO)  .BZINdOO)  ,DCdOO)  ,DGdOO)  ,DDd00)  ,Pd00)  ,RHOdOO)  , 
&ANUdOO)  ,CAdOO)  ,QTCAdOO)  ,OMG(300)  ,RTPdOO)  ,GP  (100)  ,BETAdOO)  , 
SAZOTdOO)  ,BZOT(100)  ,REY(100)  ,Q(100)  , BET (100)  ,GR(300)  ,GI(300)  , 
&DB(300) ,OMGX(300) ,OMGP(300) ,PHASE(300) ,PH(300) 

INTEGER  TURB,CASE 
CHARACTER  TITLE*50 ,THEOR*50 
DATA  PI, TPI/3. 1415926, 6. 2831853/ 

DATA  AMU, RE, GAM, SIG/0.2679E-08, 247104. ,1.4,0.8414/ 

DATA  N/10/ 

C 

C  TO  NEGATE  THE  PLOT  ROUTINE  SET  ICAS  TO  ZERO. 

C 

Qitit  it  -kit  it  If  *  -k  if  if  ft  it  it  *  -kit  if  if  -k  it  it  If  it  it  it  it  it  it  it  if  it  it  it  it  it  it  it  ft  it*  it  it  it  it  it  *  it  It  ii  it  it -kit  it  it  it  it  it  if -k  if  it*  it  it  it*  it* 


c  * 

C  THIS  PROGRAM  IS  DESIGNED  TO  CALCULATE  THE  THEORETICAL  * 

C  GAIN  AND  PHASE  SHIFT  FOR  A  GIVEN  PNEUMATIC  TRANSMISSION  * 

C  LINE.  THE  PROGRAM  FIRST  CALCULATES  THE  PRESSURE  DROP  * 

C  ALONG  THE  LINE  BY  TREATING  THE  LINE  AS  A  SERIES  * 

C  OF  SHORTER  SEGMENTS.  WITH  THE  INPUT  PRESSURE  FOR  * 

C  EACH  SEGMENT,  THE  PROPAGATION  OPERATOR  IS  CALCULATED  * 

C  USING  THE  LAMINAR  THEORY  OR  A  SELECTED  TURBULENT  * 

C  THEORY.  THE  SELECTION  OF  THE  THEORY  TO  BE  USED  * 

C  AND  THE  SELECTION  OF  THE  ORIFICE  CONDITIONS  ARE  * 

C  INTERACTIVE.  THE  PROGRAM  PROVIDES  A  TABULAR  * 

C  OUTPUT  ON  FILE  'ACROUT'  WHICH  IS  DESCRIBED  BELOW.  * 

C  THE  PROGRAM  ALSO  PRODUCES  TWO  PLOTS:  THEORETICAL  * 

C  AND  EXPERIMENTAL  GAIN  VS.  FREQUENCY,  AND  * 

C  THEORETICAL  AND  EXPERIMENTAL  PHASE  SHIFT  VS.  * 

C  FREQUENCY.  * 

C  * 

Q*********************************************************************** 
Q* **************************************** ****************************** 

C  * 

C  ACROUT  IS  THE  TABULAR  OUTPUT  FILE  OF  THE  FOLLOWING  * 

C  * 

C  PRESENTED  BY  LINE  SEGMENT  * 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

Q**** 

c 

c 

c 

c 

c 

c 

c 


Q**** 

c 

c 

c 

c 

c 

c 

c 

c 

c 


LENGTH 

DIAMETER 

PRESSURE 

ADIABATIC  SPEED  OF  SOUND 
FUNDAMENTAL  FREQUENCY 
DENSITY 

REYNOLDS  NUMBER 

THEORETICAL  RESULTS 
FREQUENCY 

NON-DIMENSIONAL  FREQUENCY 
NON-DIMENSIONAL  ATTENUATION 
PHASE  SHIFT 
GAIN 

TURBULENCE  INDICATOR 

EXPERIMENTAL  RESULTS 

FREQUENCY 
SENDING  PRESSURE 
RECEIVING  PRESSURE 
GAIN 

PHASE  SHIFT 

*********************************** A*************************** 

ATTEN  IS  AN  OUTPUT  FILE  OF  THE  NON-DIMENSIONAL 
FREQUENCY  AND  ATTENUATION  USED  AS  AN  INPUT 
FILE  FOR  THE  PLOT  ROUTINE  'PLOTM*  TO  CREATE 
PLOTS  AS  PRODUCED  BY  BROWN  OF  NON-DIMENSIONAL 
FREQUENCY  VS  NON-DIMENSIONAL  ATTENUATION 

OPEN (11 ,FILE=' D20000 ' ) 

REWIND  11 

OPEN (12 , F ILE= ' ACROUT ' ) 

REWIND  12 

OPEN  (13, FILE=' ATTEN ' ) 

REWIND  13 

OPEN ( 14, FILE=' GAIN ' ) 

REWIND  14 

CALL  PLOTS (0. ,0. ,9) 

★★★★**★★★*★**★**★***★★★****★***★★★*★**★*★★★***********★****★*★< 

TAPE  11  IS  THE  INPUT  FILE  CONTAINING  THE  FOLLOWING: 

ICAS  IS  THE  CASE  NUMBER 

OP  IS  THE  LINE  CONDITION;  1.0  FOR  OPEN  LINE, ELSE  0.0 
TR  IS  THE  AMBIENT  TEMPERATURE  IN  DEGREES  RANKINE 
PG  IS  THE  BAROMETRIC  PRESSURE  (PSIA) 

PT  IS  THE  STATIC  LINE  PRESSURE  (PSIG) 

FL  IS  THE  MASS  FLOW  RATE  (LBM/SEC) 

D1  IS  THE  OUTSIDE  LINE  DIAMETER  (IN) 

DORF2  IS  THE  SQUARE  OF  THE  ORIFICE  DIAMETER  (IN* *2) 


C  DI  IS  THE  LINE  SEGMENT  DIAMETER  (IN)  * 

C  AD  IS  THE  LINE  SEGMENT  LENGTH  (IN)  * 

C  LINE  SEGMENTS  1  &  2  ARE  THE  DOWNSTREAM  DYNAMIC  * 

C  PRESSURE  TRANSDUCER  BRANCH.  * 

C  THE  UPSTREAM  END  OF  THE  TEST  SECTION  IS  DENOTED  * 

C  AS  SEGMENT  NUMBER  10.  * 

C  * 

c  * 

C  THE  FOLLOWING  VARIABLES  ARE  USED  IN  THIS  PROGRAM  * 

C  CT  IS  THE  RATIO  OF  TURBULENT /ADIABATIC  CAPACITANCE  * 

C  TL  IS  THE  INERTANCE  CONSTANT  * 

C  ZEND  IS  THE  REAL  PART  OF  THE  TERMINAL  IMPEDANCE  (PSI/CIS)  * 

C  ZENDI  IS  THE  IMAGINARY  PART  OF  THE  TERMINAL  IMPEDANCE  (PSI/CIS)* 

C  AMU  IS  THE  DYNAMIC  VICOSITY  (PSI-SEC)  * 

C  RE  IS  THE  GAS  CONSTANT  FOR  AIR  (IN/SEC-R)  * 

C  GAM  IS  THE  RATIO  OF  SPECIFIC  HEATS  * 

C  SIG  IS  THE  SQUARE  ROOT  OF  THE  PRANDTL  NUMBER  * 

C  N  IS  THE  NUMBER  OF  LINE  SEGMENTS  * 

C  REYT  IS  THE  REYNOLDS  NUMBER  * 

C  * 

^******]lr*****«***ilt*********ilr  *******  ************************************* 

READ(ll,*)ICAS,OP 
READ(11,*)TR,PG,PT,FL,D1 
IP(OP.EQ.l. )GOTO  110 
READ(ll,*)DORF2 
110  CONTINUE 

DO  898  1=1, N 
READ(11,*)DI(I) ,AD(I) 

898  CONTINUE 

IF(OP.EQ.l.)GOTO  100 


C  IF  DORF2  =0.  THEN  THE  DATA  IS  FOR  A  BLOCKED  LINE  CASE.  * 

C  IN  THIS  CASE  THE  MODEL  SELECTION  IS  BYPASSED  AND  THE  * 

C  CALCULATIONS  ARE  DONE  USING  THE  LAMINAR  THEORY.  * 

C  CASE  IS  SET  TO  3  * 

C  N1  IS  SET  TO  3  * 

C  Nl=3  SELECTS  THE  CORRECT  HEADING  FOR  THE  OUTPUT  PILE.  * 

C  * 

Qit'k’kititicitliititililitilicieitililicititifit'k'kil’kifkititieitit'kifkitititititifitititicitit'ititicifkitit'kififk'kirltititicitititit 

IF(DORF2.EQ.0. )  THEN 
CASE=3 
Nl=3 
GOTO  96 
ENDIF 

100  CONTINUE 
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Q********************* ******************************************** ****** 

c  * 

C  THIS  INTERACTIVE  SECTION  ALLOWS  THE  SELECTION  OF  * 

C  A  LAMINAR  OR  TURBULENT  MODEL.  IT  ALSO  ALLOWS  THE  * 

C  SELECTION  OF  ONE  OF  TWO  LOAD  IMPEDANCE  MODELS.  * 

C  CASE  IS  THE  MODEL  INDICATOR  FOR  THE  LAMINAR  OR  * 

C  TURBULENT  MODEL.  IF  A  COMBINATION  OF  LAMINAR  AND  * 

C  TURBULENT  MODELS  IS  SELECTED,  THE  PROGRAM  WILL  * 

C  SWITCH  TO  THE  TURBULENT  MODEL  IF  THE  REYNOLDS  NUMBER  * 

C  IS  GREATER  THAN  2300  AND  THE  TURBULENT  ATTENUATION  * 

C  IS  GREATER  THAN  THE  LAMINAR  VALUE.  OTHERWISE  THE  LAMINAR  * 

C  MODEL  IS  USED.  * 

C  N1  IS  THE  LOAD  IMPEDANCE  INDICATOR.  YOU  HAVE  A  CHOICE  * 

C  OF  EITHER  CHOKED  OR  NON-CHOKED  FLOW  THROUGH  THE  * 

C  ORIFICE.  * 

C  * 

c 

PRINT* , ’  ' 

PRINT*,"  ' 

PRINT*,' 

PRINT* , '  ' 

PRINT* , '  ' 

PRINT*,'  PLEASE  ENTER  VALUE  OF  CASE  ' 

PRINT* , '  ' 

PRINT*,'  ' 

PRINT* , '  CASE ' 

PRINT* , '  ' 

PRINT*,'  1  LAMINAR/MOD.  CONST.  LRC 

PRINT* , '  2  LAMINAR /MODIFIED  LRCG ' 

PRINT* , '  3  LAMINAR  THEORY  ONLY ' 

PRINT*,'  ' 

PRINT*,'  ' 

PRINT* , ' 

PRINT*,' 

PRINT*,' 

PRINT*, '  ' 

PRINT* , '  ' 

PRINT*,'  ' 

PRINT*,'  ' 

PRINT* , ' ?  CASE ' 

PRINT* , '  ' 

READ*, CASE 

IF(OP.EQ.l.O)  GOTO  96 
PRINT*,'  ' 

PRINT*,'  ' 

PRINT*,'  ' 

PRINT*,'  PLEASE  ENTER  VALUE  OF  N1 ' 

PRINT*, '  ' 


o  n 


PRINT*, '  Nl' 

PRINT* , '  ' 

PRINT*,'  1  CHOKED  FLOW  IN  THE  ORIFICE' 

PRINT*,'  2  NON-CHOKED  PLOW  IN  THE  ORIFICE' 

PRINT* , '  ' 

PRINT*, '  ' 

PRINT*,'?  Nl' 

PRINT* , '  ' 

READ*,N1 

96  WRITE(12,612)ICAS 

612  FORMAT ( IX, TIO, 'RESULTS  OF  RUN  NUMBER’ ,T32 , 17) 

IF ( CASE. EQ. 3) THEN 
WRITE (12, 804) 

804  FORMAT (IHO, TIO, 'CALCULATED  USING  LAMINAR  THEORY  THROUGHOUT’) 

ELSE 

IF ( CASE. EQ.l) THEN 
WRITE (12, 80 7) 

807  FORMAT (IHO, TIO, 'CALCULATED  USING  MOD.  CONST.  LRC  THEORY  IN  THE  TUR 
&BULENT  REGION') 

ELSE 

WRITE (12, 808) 

808  FORMAT (IHO, TIO, 'CALCULATED  USING  THE  MODIFIED  LRCG  THEORY  IN  THE 
&TURBULENT  REGION') 

END  IF 

IF(OP.EQ.l.O)  GOTO  500 
1F(N1.EQ.  DTHEN 
WRITE (12, 809) 

809  FORMAT (lx, TIO, 'REFLECTION  COEFFICIENT  CALCULATED  ASSUMING  CHOKED 
&FLOW  IN  THE  ORIFICE') 

ELSE 

WRITE (12, 820) 

820  FORMAT (lx, TIO, 'REFLECTION  COEFFICIENT  CALCULATED  ASSUMING  NON- 
&CHOKED  FLOW  IN  THE  ORIFICE’) 

END  IF 

500  CONTINUE 
END  IF 
C 

THE  FOLLOWING  VARIABLES  ARE  USED  AS  PROGRAM  COUNTERS. 

ICT=1 
IND=1 
NNN=1 
C 

C  OP=0  SIGNIFIES  BLOCKED  OR  RESTRICTED  FLOW. 

C  CALCULATES  INITIAL  REYNOLDS  NUMBER  BASED  ON  THE  LINE  DIAMETER 
C  UPSTREAM  OF  THE  TEST  SECTION. 

C  (4/(AMU*PI*32.2*12) )  =  1.22999E06 

C 

IF(OP.EQ.O.)  THEN 
REYT=1. 22999E+06*FL/D1 
ELSE 
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C  TOTL  IS  THE  LINE  LENGTH  (INCHES)  FROM  THE  STATIC 
C  PRESSURE  TRANSDUCER  TO  THE  TERMINATION  OF  THE 

C  LINE.  THE  DISTANCE  FROM  THE  STATIC  TRANSDUCER 

C  TO  THE  BEGINNING  OF  THE  TEST  SECTION  IS  2.0  INCHES. 

TOTL=2 . 0 
DO  901  1=3, N 
TOTL=TOTL+AD(I) 

901  CONTINUE 

REYT=PT* (PT/2 . +PG) * (DI(N) **3 . ) / (32 . * (AMU**2 ) *TOTL*RE*TR) 
IF(REyT.LT.2300.)  GOTO  902 
REVT=(REyT*32./.1582)** (4./7.) 

902  CONTINUE 
END  IF 
PT=PT+PG 

C 

Q*********************************************************************** 

C  * 

C  THE  PRESSURE  DROP  BETWEEN  THE  STATIC  PRESSURE  * 

C  TRANSDUCER  AND  THE  BEGINNING  OF  THE  TEST  LINE  IS  * 

C  CALCULATED.  P (N)  IS  THE  LINE  PRESSURE  AT  THE  UPSTREAM  * 

C  END  OF  THE  TEST  LINE.  PF  IS  THEN  SET  EQUAL  TO  P  (N) .  * 

C  * 

IF(REyT.LT.2300. )  THEN 

DP=PT-SQRT (PT**2 . +2 . 27E-10*REyT*TR/ (Dl**3. ) ) 

ELSE 

DP=PT-SQRT (PT**2 .+1.12226E-12*REyT**l. 75*TR/ (Dl**3. ) ) 

END  IF 

P(N)=PT+DP-PG 

PF=P(N) 

IP(REyT.EQ.O. )GOTO  999 
IF(REyT.LT.2300)THEN 
F=64. /REYT 
ELSE 

F=0 . 3164/ (REyT**0 .25) 

END  IF 
F2=SQRT (F) 

TL=((F2+l.)*((F2+2.)**2))/(4.*((2*F2)+1.0)) 

CT=GAM- ( (GAM-1 . ) /TL) 

C  THIS  PROGRAM  USES  THE  TERM  (CT*AGM(I) )  * 

C  THIS  IS  THE  TURBULENT  CAPACITANCE.  IF  THE  * 

C  ADIABATIC  CAPACITANCE  IS  TO  BE  USED,  SET  CT=1.0  * 

C  AFTER  THIS  POINT.  * 

C  THE  TERM  (TL*RHO ( I ) /AR (I ) )  IS  THE  TURBULENT  * 

C  INERTANCE.  FOR  ADIABATIC  INERTANCE  SET  TL=1.0  * 

999  CONTINUE 

DO  899  1=1, N-2 
IF(OP.EQ.l. )  GOTO  895 


134 


c 

C  CALCUIATES  THE  REYNOLDS  NUMBER  IN  THE  TEST  SECTION. 

C 

REYT=1.22999E+06*FL/DI(10) 

895  CONTINUE 
C 

Q'*  *****'**  ♦***  ir  ****  **  *****!*  4:  *  Ik  **  ik  ******  sir  **  it  *  ilr  *******  jkr  ******  4r  ******  * 

c  * 

C  THE  PRESSURE  DROP  ACROSS  EACH  SEGMENT  OF  THE  LINE  IS  * 

C  NOW  CALCULATED.  THE  INDEX  I  IS  SET  FROM  1  TO  N- 2  TO  * 

C  SKIP  THE  DOWNSTREAM  TRANSDUCER  BRANCH.  THE  INDEX  K  * 

C  IS  USED  TO  REVERSE  THE  ORDER  OP  THE  SEGMENTS  SINCE  #10  * 

C  IS  THE  UPSTREAM  END  OF  THE  LINE.  * 

C  * 

^*********************************************************************** 
C  PF=PF+PG 

K=N+1-I 

IF(REYT.LT.2300.)  THEN 

DO=PF-SQRT (PF**2 . +1. 135E-10*AD(K) *REYT*TR/ (DI (K) **3. ) ) 

ELSE 

DP=PF-SQRT (PF**2 .+5 .61128E-13*REYT**1. 75*AD (K) *TR/ 

&  (DI(K)**3.)) 

ENDIF 

P(N-I)=PF+DP-PG 

L=(N-I) 

WRITE  (*,  66  )L,P(L) 

66  FORMAT  (T5, 'PC  ,11,')=  ',F10.7) 

PF=P (N-I) 

C 

C  PEND  IS  THE  LINE  PRESSURE  AT  THE  DOWNSTREAM  END  OF  THE 

C  TEST  SECTION. 

C 

PEND=P(2) 

IF(OP.EQ.l.)  PEND=0.0 
899  CONTINUE 

C  THE  PRESSURE  IN  THE  DOWNSTREAM  DYNAMIC  TRANSDUCER  BRANCH 

C  IS  SET  EQUAL  TO  THE  PRESSURE  IN  THE  LINE  AT  THE  BRANCH  POINT. 

P(2)=P(5) 

P{1)=P(2) 

C  DETERMINE  REYNOLDS  NUMBER  FOR  A  LINE 

TOTL=0 . 0 
DO  2  1=3, N 
TOTL=TOTL+AD(I) 

2  CONTINUE 

C  RHO(N)  IS  THE  DENSITY  AT  THE  BEGINNING  OF  THE  TEST  LINE. 

RHO(N)=(PG+P (N) ) / (RE*TR) 

REYT= (P (N) -PEND) * { (P (N) +PEND) /2 . +PG) * (DI(N) **3 . ) /32 . /AMU/AMU 
REYT=REYT /TOTL/ RE /TR 
IF  (REYT.LT.2300.)  GOTO  3 
REYT= (REYT*32 . /. 1582 ) ** (4 . /7 . ) 


on  on  ooooooonoonoooooooo  on 


QM  IS  THE  MASS  FLOW  RATE  (LBF-S/IN) 

QM*  386 . 04=^M (LBM/S ) 

3  QM=REYT*P1*DI (N) *AMU/4 . 

QM1=386.04*QM 
WRITE (*, 67) REYT 

67  FORMAT (T5 THE  REYNOLDS  NUMBER  IS  ’,F8.2) 

WRITE (*,68)QM1 

68  FORMAT (T5,' THE  MASS  FLOW  RATE  IS  ',1PE11.4,'  LBM/SEC’) 
*********************************************************************** 


* 

* 

THE  VARIABLES  IN  THE  FOLLOWING  SECTION  ARE  DEFINED  AS;  * 

RHO=DENSITY  CALCULATED  BY  THE  IDEAL  GAS  LAW  (I£F-S**2/IN**4)  * 

ANU=KINEMATIC  VISCOSITY  (IN**2/S)  * 

CA=ADIABATIC  SPEED  OF  SOUND  (IN/S)  * 

AR=  CROSS-SECTIONAL  AREA  OF  LINE  (IN* *2)  * 

CV=VISCOUS  CHARACTERISTIC  FREQUENCY  (OMEGA-NU)  (RAD/S)  * 

OMT=THERMAL  CHARACTERISTIC  FREQUENCY  (RAD/S)  * 

CNA=VISCOUS  RESISTANCE  (LBF-S/IN**6)  * 

Q=VOLUMETRIC  FLOW  RATE  (IN**3/S)  * 

AGM=ADIABATIC  CAPACITANCE  * 

RVT=LINE  RESISTANCE  * 

FN=FUNDAMENTAL  FREQUENCY  (RAD/S)  * 

PBR=LINE  SEGMENT  PRESSURE  (PSIA)  * 

* 


*********************************************************************** 

DO  23  1=1, N 
PBR=P(I)+PG 
RHO(I)PBR/(RE*TR) 

ANU(I)=AMU/RHO(I) 

CA ( I ) =SQRT (PBR*GAM/RHO ( I ) ) 

AR(I)=PI*DI (I) **2 ./4. 

CV(I)=(8. *PI*ANU(I))/AR(I) 

OMr(I)=CV(I) /(SIG*SIG) 

CNA(I)=(8.*PI*AMU) /(AR(I)*AB(I)) 

Q(I)=QM/RHO(I) 

23  CONTINUE 

THE  FLOW  IN  THE  DYNAMIC  TRANSDUCER  BRANCH  IS  SET 
EQUAL  TO  ZERO. 

Q(1)=0.0 
Q(2)=0.0 
GMNI=.5* (GAM-1. ) 

DO  26  1=1, N 

REY(I)=(4.*RHO(I)*Q(I) )/(PI*DI (I)*AMU) 

PBR=P(I)+PG 
AGM(I)=AR(I) / (GAM*PBR) 

QTCA(I)=0.25*CA(I) 

RVT(I)  IS  THE  DC  LINE  RESISTANCE 


nnnnonoonoooooononnoooooonooo  oo 


1.75*RVT(I)  IS  THE  AC  (INCREMENTAL)  LINE  RESISTANCE 

RVT(I)=.3164*AMU* (REY (I) ** .75) /(2 . *AR (1) * (DI (I) **2. ) ) 

RVT(I)=RVT(I)*1.75 

PN(I)=QTCA(I)/AD(I) 

26  CONTINUE 
NST=1 
DW=0. 

Y=l. 

M=0 


**** ********************************************  ****************  ******* 

LOOP  80  CALCULATES  THE  THEORETICAL  VALUES  FOR  THE  * 

PLOT  ROUTINE.  THE  MODIFIED  EQUATIONS  OF  KRISHNAIYER  * 

AND  LECHNER  ARE  USED  FOR  BLOCKED  AND  LAMINAR  FLOW.  * 

THE  PROGRAM  ALSO  CALCULATES  THE  PROPAGATION  CONSTANT  * 

BASED  ON  THE  SELECTED  TURBULENT  MODEL  AND  COMPARES  * 

THE  LAMINAR  AND  TURBULENT  VALUES.  IT  SELECTS  THE  * 

GREATER  VALUE  FOR  USE  IN  THE  OUTPUT  ROUTINE.  * 

* 

* 

THE  VARIABLES  USED  ARE:  * 

* 


RN=REAL  PART  OF  Z  (EQN  6  WILKIN)  * 
DC=COEFFICIENT  (EQN  11  WILKIN)  * 
DGCOEFFICIENT  (EQN  10  WILKIN)  * 
DD=DENOMINATOR  OF  EQN  7  (WILKIN)  * 
GN=REAL  PART  OF  Y  (EQN  7  WILKIN)  * 
ALN=IMAGINARY  PART  OF  Z  (EQN  6  WILKIN)  * 
CN=I MAG INARY  PART  OF  Y  (EQN  7  WILKIN)  * 
TEM1=REAL  PART  OF  ZY  * 
TEM2=IMAGINARY  PART  OF  ZY  * 
AN=REAL  PART  OF  SQRT(ZY)  * 
BTN= IMAGINARY  PART  OF  SQRT(ZY)  * 
OX=NONDIMENSIONAL  FREQUENCY  * 
OY=NONDIMENSIONAL  ATTENUATION  * 
TURB=TURBULENT  MODEL  INDICATOR  * 


Q*** ******************************************  ************************** 

40  DO  80  J=1,NST 
M=M+1 
C 

C  INCREMENT  THE  FREQUENCY  BY  5  HZ  AND  CONVERT  TO  RAD/SEC. 

C 

Y=Y+DW 
W=TPI*Y 
DO  27  1=1, N 
ARG=.5*SQRT(W/CV(I) ) 

RN(I)=CNA(I) *(.375+ARG+(. 375/(4. *ARG) ) ) 


DC (I) = . 25+SQRT (W/OMT {!))+. 125*SQRT (OMT ( I) /V) 

DG (I) =SQRT (W/OMT (I) ) - . 125*SQRT (OMT (I) /W) 

DD (I)  =DC (I)  *DC (I)  +DG (I)  *DG (I) 

GN ( I ) =W* ( GAM- 1 . ) *AGM ( I ) * DG ( I ) /DD ( I ) 

ARG= . 5**SQRT (CV ( I) /W) 

ALN ( 1 ) =RHO ( I ) * ( 1 . +ARG- ( ARG* ( 1 5 . *C V ( I ) / ( W* 6 4 . ) ) ) ) /AR ( I ) 

TURB=0 

CN ( I) =AGM(I) * (1.+ ( (GAM-1 . ) *DC (I) /DD (I) ) ) 

TEMP=-W*W 
LT=TL*RHO { I ) /AR ( I ) 

DENOM=RVT (I) **2+ ( (W**2) * (LT**2) ) 

Q*** ************************************************** ****************** 

C  THIS  SECTION  CALCULATES  THE  PKDPAGATION  * 

C  OPERATOR  SQRT(ZY).  THE  THEORY  SELECTED  AT  * 

C  THE  BEGINNING  OF  THE  PROGRAM  IS  USED  IN  * 

C  THIS  SECTION.  * 

C  410  MODIFIED  CONSTANT  LRC  * 

C  420  MODIFIED  LRCG  * 

C  430  LAMINAR  * 

Q********************************* ******* ******************************* 

GOTO ( 410 , 420 , 430 ) CASE 
410  TEM1=RN  (I)  *GN  (I)  +TEMP*ALN  (I)  *CN  (I) 

TEM2=W* (RN (I) *CN (I) +GN (I) *ALN (I) ) 

CALL  RTCMP(ARG1,ARG2,TEM1,TEM2) 

AN(I)=ARG1 

BTN(I)=ARG2 

TEM1=TEMP'*TL* RHO  ( I)  /AR  ( I)  *CT»AGM ( I ) 

TEM2  =W*RVT ( 1 ) *CT  * AGM ( I ) 

CALL  RTCMP(ARG1,ARG2,TEM1,TEM2) 

TAN=ARG1 

TBTN=ARG2 

IF(REY(I) .LE.2300)  GOTO  27 
IF (AN (I) .GE. TAN)  GOTO  27 
AN(I)=TAN 
BTN(I)=TBTN 
TURB=1 
GOTO  27 

420  TEM1=RN (I) *GN (I) *TEMP*ALN (I) *CN (I) 

TEM2=W* (BN(I)*CN(I) *GN (I) *ALN (I) ) 

CALL  RTCMP(ARG1,ARG2,TEM1,TEM2) 

AN(I)=ARG1 

BTN(I)=ARG2 

CT=GAM- ( (W**2) * (GAM-1 . ) *LT*RHO (I) / (AR ( 1) *DENOM) ) 

TEM1=( ( (RVT(I) *(GAM-1.) *AGM(I)*RHO(I) * (W**2) ) /(AR(I) *DENOM) ) 

&) )+(TEMP*TL*RHO(I) *CT*AGM(I) /AR(I) ) 

TEM2=( (RVT(I)* (GAM-1.) *AGM(I)*RHO(I) *LT* (W**3) ) / (AR(I) *DENOM) ) 
&+(W*RVT(I)*CT*AGM(I) ) 

CALL  RTCMP(ARG1,ARG2,TEM1,TEM2) 

TAN=ARG1 

TBTN=ARG2 

IF (REY( I) .LE.2300)  GOTO  27 
IF  (AN (I)  .GE. TAN)  GOTO  27 


A1J(I)=TAN 
BTN(I)=TBTN 
TURB=1 
GOTO  27 

430  TEM1=RN(I)*GN(I)+TEMP*ALN(I)*CN(I) 

TEM2  =W*  (RN  ( I )  *CN  ( I ) +GN  ( I )  *AIiN  ( I )  ) 

CALL  RTCMP(ARG1,ARG2,TEM1,TEM2) 
AN(I)=ARG1 
BTN(I)=ARG2 
27  CONTINUE 

OX=AR ( 10 ) *W/ANU ( 10 ) /PI 
OY=AR ( 10 ) *CA ( 10 ) * AN ( 10 ) /PI /ANU ( 10 ) 
WRITE(13,75)OX,OY 
75  FORMAT(T5,F10.5,T20,F10.5) 

C  CALCULATE  ZO 
DO  32  1=1, N 
TEM1=W*ALN(I) 

TEM2=W*CN(I) 

TEM3=RN(I) 

TEM4=GN(I) 

CALL  CMPDV(ARG1,ARG2,TEM3,TEM1,TEM4,TEM2) 
CALL  RTCMP(AZRNI,BZRNI,ARG1,ARG2) 

AZRN (I) =AZRN1 
BZRN(I)=BZRNI 
IF{TURB.EQ.O)  GOTO  32 
TEM1=RH0(1) * (P (I) +PG) /AR (I) /AR (I) 
TEM2=-RVT(T) *(P(I)+PG) /W/AR(I) 

CALL  RTCMP(ARG1,ARG2,TEM1,TEM2) 
A2RN(I)=ARG1 
BZRN(I)=ARG2 
32  CONTINUE 
C  CALCULATE  Z  IN  1 
1=1 

TEMP=AN(I)*AD(I) 

IF (TEMP, GT. 88.)  GO  TO  80 
ARG1=C0SH (TEMP) 

ARG2=SINH(TEMP) 

TEM5=BTN(I)*AD(I) 

TEMP=COS(TEM5) 

TEM1=ARG1*TEMP 

TEM3=ARG2*TEMP 

TEMP=SIN(TEM5) 

TEM2=ARG2*TEMP 

TEM4=ARG1*TEMP 

CALL  CMPDV ( ARG 1 , ARG2 , TEMl , TEM2 , TEM3 , TEM4 ) 
TEM1=AZRN(I) 

TEM2=BZRN(I) 

CALL  CMPMP ( TEM3 , TEM4 , TEMl , TEM2 , ARGl , ARG2 

AZIN(I)=TEM3 

BZIN(I)=TEM4 


C  CALCULATE  Z  IN  3 

C  SEGMENT  3  IS  THE  DOWNSTREAM  END  OP  THE  TEST  SECTION 


1=3 

TEM1=AZRN(I) 

TEM2=BZRN(I) 

C 

^  **  A  *  4k  *  ik  ******  *4r  ir  Ik  ***  **************************************  ********  * 


C  CALCULATES  THE  TERMINAL  IMPEDANCE:  IF  OP=l  (UNRESTRICTED 

C  FLOW) ,  SETS  REAL  PART  EQUAL  TO  ZERO 

C  CALCULATES  IMAGINARY  PART 

C  IF  OP=0  (RESTRICTED  FLOW)  ,  SETS  IMAGINARY  PART  TO  ZERO 

C  CALCULATES  REAL  PART 

C  FOR  BLOCKED  LINE  RE=  INFINITY 

C  FOR  ORIFICE  FLOW  RE  CALCULATED  USING  REFLECTION 

C  COEFFICIENTS 

C 


* 

* 

* 

* 

* 

* 

* 


* 


Q*********************************************************************** 


C  IF  (OP.EQ.l.)  THEN 

ZEND=0.0 

ZENDI=-W*.4*RHO(I) *DI(I) /AR(I) 

ELSE 

ZENDI=0.0 

IF (DORF2 . EQ . 0 . ) THEN 

ZEND=1.0E99 

ELSE 

SQD3=DI(I) **2 

C  Nl  IS  AN  INDICATOR  TO  SELECT  THE  METHOD  OF  CALCULATING 

C  THE  REFLECTION  COEFFICIENT. 

C  Nl=l  CHOKED  FLOW  IN  THE  ORIFICE 

C  Nl=2  NON-CHOKED  FLOW  IN  THE  ORIFICE 
C 

IF(N1.EQ.1.)THEN 
RC= (SQD3-DORF2 ) / (SQD3+DORF2) 

ELSE 

RC= (SQD3-2 . *DORF2) / (SQD3+DORF2) 

END  IF 

RCC=(1.+RC) /(l.-RC) 

ZEND=4 . *RHO ( I ) *CA ( I ) *RCC/ ( PI*SQD3) 

END  IF 
ENDIF 
TEM3=ZEND 
TEM4=ZENDI 
TEM5=AN(I) 

TEM6=AD(I) 

TEM7=BTN(I) 

CALL  CALZ IN { AARG , BARG , TEMl , TEM2 , TEM3 , TEM4 , TEM  5 , TEM6 , TEM7 ) 

AZIN(I)=AARG 

BZIN{I)=BARG 

C  CALCULATE  Z  IN  4  AND  5 
DO  34  1=4,5 
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TEM1=A2RN(I) 

TEM2=BZRN(I) 

TEM3=AZIN{I-1) 

TEM4=BZIN(I-1) 

TEM5=AN{I) 

TEM6=AD(I) 

TEM7=B1N (I) 

CALL  CALZ IN ( AARG , BARG , TEMl , TEM2 , TEM3 , TEM4 , TEM5 , TEM6 , TEM7 ) 
AZIN(I)=AARG 
34  BZIN(I)=BARG 
C  CALCULATE  Z  IN  2 
1=2 

TEM1=AZRN(I) 

TEM2=BZRN(I) 

TEM3=AZIN(I-1) 

TEM4=BZIN (I-l) 

TEM5=AN(I) 

TEM6=AD(I) 

TEM7=BTN(I) 

CALL  CALZIN (AARG, BARG , TEMl ,TEM2 ,TEM3 , TEM4 ,TEM5 ,TEM6 ,TEM7 ) 

AZIN(I)=AARG 

BZIN(I)=BARG 

C  CALCULATE  RECEIVING  Z  FOR  LINE  6 
1=^ 

TEMl =AZIN (I-l) 

TEM2=BZIN(I-1) 

TEM3=AZIN(I-4) 

TEM4=BZIN(I-4) 

CALL  ZEBRA (AZOTI,  BZOTI,TEMl,TEM2,TEM3,TEM4) 

AZOT(I)=AZOTI 

BZOT(I)=BZOTI 

C  CALCULATE  Z  IN  6  INCLUDING  R-TRANSDUCER 
TEM1=AZRN(I) 

TEM2=BZPN(I) 

TEM3=AN(I) 

TEM4=AD(I) 

TEM5=BTN(I) 

CALL  CALZIN (AZINI , BZINl , TEMl , TEM2 , AZOTI , BZOTI , TEM3 , TEM4 , TEM5 ) 

AZIN(I)=AZINI 

BZIN(I)=BZINI 

C  CALCULATE  Z  IN  7,8,9  AND  10  INCLUDING  R-TRANSDUCERS 
DO  39  1=7,10 
TEM1=AZRN(I) 

TEM2=BZRN(I) 

TEM3=AZIN(I-1) 

TEM4=BZIN(I-1) 

TEM5=AN(I) 

TEM6=AD(I) 

TEM7=BTN(I) 

CALL  CALZIN (AZINI, BZINl, TEMl ,TEM2 ,TEM3 ,TEM4 ,TEM5 ,TEM6 ,TEM7) 
AZIN(I)=AZINI 


39  BZIN(I)=BZINI 

C  CALCULATE  P5/P6  INCLUDING  R-TRANSDUCER 
1=6 

TEMP=BTN(I) *AD(I) 

CSBI1=C0S (TEMP) 

SNBI1=SIN (TEMP) 

TEMP=AN(I)*AD(I) 

ARG1=C0SH (TEMP) 

ARG2=SINH(TEMP) 

TEM1=ARG1*CSBI1 

AZOTI=AZOT{I) 

BZOTI=BZOT(I) 

AZINI=AZIN(I) 

BZINI=BZIN(I) 

CALL  CMPDV(TEM7 ,TEM8,AZ0TI,BZ0TI,AZINI,BZINI) 

CALL  CMPMP(TEM3,TEM4,TEM7,TEM8,TEM1,TEM2) 

TEM5=ARG2*CSBI1 

TEM6=ARG1*SNBI1 

TEM9=AZRN(I) 

TEM10=BZRN(I) 

CALL  CMPDV(TEM7,TEM8,AZOT1,BZOTI,TEM9,TEM10) 

CALL  CMPMP(TEM1,TEM2,TEM7,TEM8,TEM5,TEM6) 

TEM1=TEM3-TEM1 
TEM2=TEM4-TEM2 
TEMP=TEM1 *TEM1+TEM2  *TEM2 
BETA ( I ) =ATAN2 (TEM2 , TEMl ) 

BET(I)=(180./PI)*BETA(I) 

RTP (I- 1) =SQRT (TEMP ) 

GP(I-1)=20.*ALOG10(RTP(I-1) ) 

C  CALCULATE  P6/P7,  P7/P8,  P8/P9,  AND  P9/P10  INCLUDING  R-TRANSDUCERS 
DO  43  1=7,10 
TEMP=BTN(I) *AD(I) 

CSBI1=C0S(TEMP) 

SNBI1=SIN(TEMP) 

TEMP=AN(I)*AD(I) 

ARGl=COSH (TEMP) 

ARG2=SINH(TEMP) 

TEM1=ARG1*CSBI1 

TEM2=ARG2*SNBI1 

AZINI1=AZIN(I-1) 

BZINI1=BZIN(I-1) 

AZINI=AZIN(I) 

BZINI=BZIN(I) 

CALL  CMPDV(TEM7,TEM8,A2INIl,BZINIl,AZINI,BZINI) 

CALL  CMPMP (TEM3 , TEM4 , TEM7 , TEM8 , TEMl , TEM2 ) 

TEM5=ARG2*CSB11 

TEM6=ARG1*SNB11 

TEM9=AZRN(1) 

TEM10=BZRN(I) 

CALL  CMPDV ( TEM7 , TEM8 , AZ INI 1 , BZ INIl , TEM9 , TEMl 0 ) 

CALL  CMPMP ( TEMl , TEM2 , TEM7 , TEM8 , TEM5 , TEM6 ) 


TEM1=TEM3-TEM1 
TEM2=TEM4-TEM2 
TEMP=TEM1*TEM1+TEM2  *  TEM2 
BETA  (I)  =ATAN2  (TEM2  ,TEM1) 

BET (I ) = (180 . /PI ) *BETA (I ) 

RTP ( I- 1 ) SQRT  (TEMP) 

43  GP(I-1)=20.*ALOG10(RTP(I-1) ) 

RrPT=RTP  (6)  *RTP  (7)  *RTP  (8)  *RTP  (9)  *RTP  (5) 

GPT=20 .  *ALCX310  (RTPT) 

BETAR=BETA (6 ) +BETA ( 7 ) +BETA ( 8 ) +BETA (9 ) +BETA ( 10 ) 

35  IF  (BETAR-LE.O. )  GOTO  36 
BETAR=BETAR-TPI 

TOGO  35 

36  IF  (BETAR.GT.-TPI)  GOTO  37 
BETAR=BETAR+TP I 

GOTO  36 

2  CONVERTS  BETAR(RADS)  TO  BETAD(DEGS) 

37  BETAD= ( 180 . /PI ) *BETAR 

2  SKIPS  WRITE  STATEMENT  FOR  OUTPUT  TABLE  HEADINGS 
IF(NNN.GT.l)  GOTO  600 
WRITE (12,810) 

810  FORMAT (1H0,T 3, 'SEGMENT' ,T16 , 'LENGTH' ,T2 9, 'DIAMETER' ,T43, 

&' PRESSURE' ,T59, 'C  ADB' ,T7 3, 'FN' ,T86, 'DENSITY' ,T98, 'REYNOLDS  NO') 

DO  803  1=1, N 
IF(OP.EQ.l.O)  GOTO  803 
IF(DORP2.EQ.0,0)  THEN 
REY(I)=0. 

END  IF 

803  WRITE (12, 802) I, AD (I) ,DI (I) ,P(I) ,CA(1) ,FN(I) ,RHO(I) ,REY(I) 

802  FORMAT(1HO,T6,I2,T11,7 (1PE12.4,2X) ) 

WRITE (12, 706) 

706  FORMAT (IH//) 

WRITE (12, 890) 

890  FORMAT (t2 5, '****************************************************') 

WRITE (12, 7 06) 

WRITE (12, 700) 

700  FORMAT (T40 , ' THEORETICAL  RESULTS ’ ) 

WRITE (12, 7 06) 

WRITE (12, 830) 

830  FORMAT ( IX, TlO, 'FREQUENCY' ,T35,' ATTENUATION ' ,T52, 'PHASE  SHIFT' ,T79 
& , ' GAIN ' , TlO 1 , ' TURBULENCE ' ) 

WRITE (12,850) 

850  FORMAT ( IX, TlOl, ' INDICATOR' ) 

WRITE (12,840) 

840  FORMAT (IX, T5 , ' (HZ) ' ,T13 , ' (NON-DIMENSIONAL) ' ,T32 , ’ (NON-DIMENSIONAL) 
&' ,T55, ' (DEG) ' ,T70, ' (DB) ' ,T82, ' (NON-DIMENSIONAL) ' ) 

WRITE (12,706) 

C  NNN  IS  THE  COUNTER  TO  SKIP  THE  HEADINGS  WRITE  STATEMENT. 

NNN=2 

600  WRITE (12, 702 )Y, OX, OY,BETAD,GPT, RTPT, TURB 


702  FORMAT (IX, T4 ,F5 . 0 ,T14 ,F10 . 5 ,T34 ,F10. 5 , T54 ,F7. 2 ,T66 , 1PE12 . 4 , T83 , IPE 
&12.4,T105,T2) 

WRITE(14,97)y,GPT 
97  PORMAT(T5,F5.0,T20,1PE12.4) 

DBT (M) =GPT 
0MB (M) =y 
PHASE (M) =BETAD 

80  CONTINUE 

GO  TO  (85,611) ,IND 
85  y=o. 

C  DW  IS  THE  FREQUENCY  INCREMENT  (HZ),  NST  IS  THE  NUMBER 
C  OF  THEORETICAL  VALUES  CALCULATED  (0-1000  HZ) 

DW=5. 

NST=200 
IND=2 
GO  TO  40 


C 

Qil^Aik**4f*******i^****«^*****«t*lti^********4r**4r*Wilr**4ri^**ilr*3lr9^****^***4r*1k******#4r 
Q* ****************************************************** **************** 
^*********************************************************************** 
^*** *************************************************** ***************** 
Q* *************************************************************** ******* 
^* ****************************************************** **************** 
^******* ******************************************************** ******** 
C  THIS  SECTION  IS  FOR  THE  CALCOMP  PLOTTER.  ONE  SET  OP  * 

C  EXPERIMENTAL  DATA  MUST  BE  INCLUDED  WITH  EACH  RUN.  * 

C  * 

C  NPTS  IS  NUMBER  OF  EXPERIMENTAL  POINTS  TO  BE  INPUT  * 

C  LSMB  IS  THE  PLOTTER  SYMBOL  TO  BE  USED  FOR  GPX  AND  PHD.  * 

C  SEE  A  CALCOMP  PLOTTER  USERS  MANUAL  FOR  DESCRIPTION.  * 

C  FREQ  IS  THE  EXPERIMENTAL  FREQUENCY  (HZ)  * 

C  PS  IS  THE  SENDING  PRESSURE,  EXPERIMENTAL  (MV)  * 

C  PR  IS  THE  RECEIVING  PRESSURE,  EXPERIMENTAL  (MV)  * 

C  PHT  IS  THE  EXPERIMENTAL  PHASE  DELAY  (MSEC)  * 

C  * 

Q*********************************************************************** 

C  NEGATES  THE  PLOT  ROUTINE  IF  ICAS  .LE.  ZERO 

611  IF (ICAS) 502, 502 ,610 
610  CALL  FACTOR(0.625) 

CALL  RECT(0. ,0. ,9.6,12.0,0. ,3) 

CALL  PLOT(.6, .6,-3) 

OMG(M+1)=0. 

OMG(M+2)=100. 

DBT(M+1) =-20. 

DBT (M+2 ) =5 . 0 

C  DRAWS  AXES  FOR  GAIN  PLOT 

CALL  AXIS(0. ,0. ,' FREQUENCY  (HERTZ) ’ ,-17,10.0,0. ,OMG(M+l) ,OMG 
&  (M+2)) 
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CALL  AXIS(0. ,0. , 'GAIN  (DECIBELS) ' ,15,8. ,90. ,DBT(M+1) ,DBT 
&  (M+2)) 

PLOTS  CURVE  FOR  THEORETICAL  GAIN 
CALL  LINE(OMG,DBT,M,1,0,4) 

READ (11, 705) TITLE 
705  FORMAT (A50) 

READ(11,*)NPTS 
WRITE (12,706) 

WRITE (12, 890) 

WRITE (12,706) 

WRITE (12,860) 

860  FORMAT (T45 EXPERIMENTAL  DATA') 

WRITE (12,207) 

207  F0RMAT(1H//,T14, 'N' ,T20, 'FREQUENCY  (HZ) ' ,T42 ,' PRESSURE  (MV) ' ,T66 
& ' GAIN  (DB) ' ,T82 , ' PHASE  SHIFT  (DEG) ' ) 

WRITE (12,870) 

870  FORMAT (1H//,T39, 'SENDING' ,T49, 'RECEIVING') 

WRITE (12, 706) 

J=0 

DO  69  I=1,NPTS 

READ ( 1 1 , * ) FREQ , PS , PR , PHT 

GPX=20 . *ALOG10 (PR/PS) 

PHD=-PHT*FREQ* . 36 
IF (PHD . LT . - 360 . ) THEN 
PHD=PHD+360. 

END  IF 

WRITE(12,209) I, FREQ, PS, PR, GPX, PHD 

209  FORMAT(1X,T12,I3,T23,F5.0,T39,F6.1,T50,F6.1,T64,F10.5,T84,F10.5) 
DB(I)=GPX 
OMGX(I)=FREQ 
IF(PHD.GT.O)  GOTO  69 
J=J+1 
PH(J)=PHD 
OMGP(J)=FREQ 
JMAX=J 
69  CONTINUE 

OMGX (NPTS+1 ) =OMG ( M+1 ) 

DB (NPTS+1) =DBT (M+1) 

OMGX (NPTS+2 ) =OMG (M+2 ) 

DB (NPTS+2 ) =DBT (M+2 ) 

PLOTS  EXPERIMENTAL  GAIN  DATA 
LSMB=0 

CALL  LINE (OMGX, DB,NPTS,1,-1,LSMB) 

CALL  SYMBOL(3. ,13. ,.165,TITLE,0. ,50) 

GOTO (111^12,1 13) CASE 

111  THEOR='MOD  CONST  LRC' 

GOTO  114 

112  THEOR*' CONST  LRCG' 

GOTO  114 

113  THEOR=' LAMINAR' 


114  CONTINUE 

CALL  SYMBOLO.  ,12 . 5,  .  165 ,THEOR,0 .  ,50) 

CALL  PLOT(15.0,-.6,-3) 

CALL  FACTOR(0.625) 

CALL  RECT(0. ,0. ,9.6,12.0,0. ,3) 

CALL  PLOT (.6,. 6, -3) 

OMG(M+1)=0. 

OMG(M+2)=100. 

PHASE (M+1) =-360. 

PHASE (M+2) =45. 

C  DRAWS  AXES  FOR  PHASE  PLOT 

CALL  AXIS (0. ,0. , 'FREQUENCY  (HERTZ) * ,-17 , 10 . ,0 . ,OMG (M+l) , 

&OMG(M+2) ) 

CALL  AXIS(0. ,0. , 'PHASE  ANGLE  (DEGREES) ' ,21,8. ,90. ,PHASE(M+1) , 
&PHASE(M+2) ) 

C  PLOTS  CURVE  FOR  THEORETICAL  PHASE 

CALL  LINE (OMG, PHASE, M, 1,0, 4) 

OMGP{JMAX+l)  =OMG(M+l) 

OMGP ( JMAX+2 ) =OMG (M+2) 

PH (JMAX+1) =PHASE (M+l) 

PH (JMAX+2) =PHASE (M+2) 

C  PLOTS  EXPERIMENTAL  PHASE  DATA 

CALL  LINE (OMGP, PH, JMAX,1,-1,LSMB) 

CALL  SYMBOLO.  ,13.  ,  .165, TITLE, 0.  ,50) 

CALL  SYMBOLO.  ,12.5,  .165, THEOR,0.  ,50) 

1003  CALL  PLOT(10.0,0,-3) 

^**********************lk**«*******1k4t**ilr*ilrilr()^***ilrikilr)li^i^4^i^ilrilri^1^ilrt^4r*ilr«lr***ltr**ilr* 

c 

501  CALL  PLOTE 

502  STOP 
END 

Qit  it  Hit  it*  it  it  It  it -kit  it  it  hit  1i  It  Ifk  It  h*  it  it  It  it  It  lilt  It  it  it  if  ic  it  It  it  it -ft  •kit'k  It  it  it  It  ic  it  it  it  It  it  it  it  It  icicit  it  it  it  it  it  ii  It  It  It  It  it  It 
Qit  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  if  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it 

C  CALCULATES  ROOT  OF  A  COMPLEX  NUMBER 
C  X  =  RE(SQRT(A+BI) 

C  Y  =  IM(SQRT(A+BI) 

C 

SUBROUTINE  RTCMP (X, Y,A,B) 

CALL  ANGL(TEMP,A,B) 

TEMP=.5*TEMP 

Y=A*A+B*B 

X=SQRT(Y) 

X=SQRT(X) 

Y=X*SlN(TEMP) 

X=X*COS (TEMP) 

RETURN 

END 


o  o  n  n  on 


c 

SUBROUTINE  CMPMP (X.Y, A1,A2 ,B1,B2) 
C  MULTIPLIES  TWO  COMPLEX  NUMBERS 
(A1+A2 J) * (B1+B2 J) =X+yj 

X=A1*B1-A2*B2 

Y=A1*B2+A2*B1 

RETURN 

END 


FINDS  THE  QUOTIENT  OF  TWO  COMPLEX  NUMBERS 
SUBROUTINE  CMPDV (Cl ,C2 ,Al ,A2 ,Bl ,B2) 

C  CALCULATES  THE  QUOTIENT  OF  TWO  COMPLEX  NUMBERS 
C  A1+A2J/B1+B2J 

TEMP=B1*B1+B2*B2 

Cl=Al*Bl+A2*B2 

C1=C1/TEMP 

C2=B1*A2-A1*B2 

C2=C2/TEMP 

RETURN 

END 

C 

Q*  ** -k  *  i(  *  it  *** -h  ***  it  **  it  * -k  it  it -k  ii  *  ^  it  ***  It  it  it  h  ■kititititifkit1t1t'k1titit1fk‘kitititit‘kit1titititiiii1iii  it  it  it  it  it 


SUBROUTINE  HS1NX(ARG,X) 

C  CALCULATES  THE  HYPERBOLIC  SINE  OF  X 
A=EXP (X) 

B=EXP (-X) 

A=A-B 
ABG=. 5*A 
RETURN 
END 
C 

Qir  it  it  -k  if  ir*  it  iHrit^t  :kir  it  it  iririr  irir  ■kit  it  ir  it  iiiHrir  flit  ik  hiHr  It  It  It  It  it  icii  it  it  it  it  it  it  it  it  it  icit  it  It  It  It  it  it  it  It  it  1r -kit  it  It  it  It  It 

c 


SUBROUTINE  HCOSX(ARG,X) 

C  CALCULATES  THE  HYBERBOLIC  COSINE  OF  X 
A=EXP (X) 

B=EXP (-X) 

A=A+B 

ARG=.5*A 

RETURN 

END 


c 


SUBROUTINE  ANGL(C,A,B) 

C  CALCULATES  THE  MAGNITUDE  OF  ANGLE  X  (IN  RADS)  AND  ASSIGNS 

C  THE  PROPER  SIGN  BASED  ON  THE  CARTESIAN  COORDINATES  (A,B) 

DATA  PI/3.1415926/ 

C=ABS (B/A) 

C=ATAN(C) 

IF  (A.GT.O.)  GO  TO  5 
IA=1 
GO  TO  7 
5  IA=0 

7  IF  (B.GT.O.)  GO  TO  10 
IB=2 

GO  TO  15 
10  IB=0 
15  IA=IA+IB+1 

GO  TO  (35,30,25,20) ,IA 
20  C=C-PI 
GO  TO  35 
25  C=-C 

GO  TO  35 
30  C=PI-C 
35  RETURN 
END 
C 

Q***** **********  ******************************************************** 


SUBROUTINE  CALZIN (AZINl ,BZINl ,AZRNl ,BZRN1 , AZIN2 ,BZIN2 , ANl ,DI1 , 
&BTN1) 

TEMP=AN1*DI1 

CALL  HC0SX(ARG1,TEMP) 

CALL  HSINX(ARG2,TEMP) 

TEMP=BTN1*DI1 

CSBI1=C0S(TEMP) 

SNBI1=SIN(TEMP) 

ZR=0. 

CALL  CMPMP ( TEMl , TEM2 , AZ IN2 , BZIN2 , ARGl , ZR ) 

CALL  CMPMP ( TEM  3 , TEM4 , AZRN 1 , BZRN 1 , ARG  2 , ZR ) 

CALL  CMPMP (TEM5 , TEM6 , AZ IN2 , BZ IN2 ,ARG2 ,ZR) 

CALL  CMPMP (TEM7 ,TEM8 , AZRNl ,BZRN1 ,ARG1 , ZR) 

Al=TEMl+TEM3 

Bl=TEM2+TEM4 

A2=TEM5+TEM7 

B2=TEM6+TEM8 

CALL  CMPMP (TEMl , TEM2 , Al , Bl , CSBIl , ZR) 

CALL  CMPMP (TEM5 ,TEM6 , A2 ,B2 , CSBIl, ZR) 

CALL  CMPMP (TEM7 , TEM8 , Al , Bl , ZR , SNBI 1 ) 

CALL  CMPMP ( TEM3 , TEM4 ,A2,B2,ZR,SNBI1) 
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TEM3=TEM3+TEM1 

TEM4=TEM4+TEM2 

TEM7=TEM7+TEM5 

TEM8=TEM8+TEM6 

CALL  CMPDV{TEM1,TEM2,TEM3,TEM4,TEM7,TEM8) 

CALL  CMPMP(AZIN1,BZIN1,TEM1,TEM2,AZRN1,BZRN1) 

RETURN 

END 

C 

(;;*********************************************************************** 

C 

SUBROUTINE  ZEBRA(C1,C2 , A1 ,A2 ,B1,B2) 

C  COMPUTES  THE  INVERSE  OF  THE  SUM  OF  THE  INVERSES  OF  TWO 
C  COMPLEX  NUMBERS . 

C  Cl+C2J=l/( (1/A1+A2J)+(1/b1+B2J) ) 

Dl=l. 

D2=0. 

CALL  CMPDV(ARG1,ARG2,D1,D2,A1,A2) 

CALL  CMPDV(ARG3,ARG4,D1,D2,B1,B2) 

ARG1=ARG1+ARG3 

ARG2=ARG2+ARG4 

CALL  CMPDV(C1,C2,D1,D2,ARG1,ARG2) 

RETURN 


Appendix  C :  Equii 

3ment 

•• 

Item 

Manufacturer 

Model 

Serial  No. 

Pneumatic  driver 

Bendix 

PC-DCC-PS 

001 

, _  i 

*'•  /* 

Amplifier/ 

ENY  A  293 

Power  Supply 

Bendix/Lambda 

-297 

Wave  Analyzers 

Hewlett  Packard 

302A 

018-02012 

018-01985 

.  I 

V  •  »  •  ♦ 

Frequency 

Computer 

726C 

ENY  0021 

Counters 

Measurement 

ENY  0035 

"  -  •  H  •' 

Dynamic  Pressure 

Kristal 

701A 

65079 

Ir'-'j-yj 

1 

Transducers 

7588 

Charge 

Kistler 

568 

704 

Amplifiers 

1429 

Digital 

Nicolet 

2090-3C 

85A06578 

1 

Oscilloscope 

_ \ 

Static  Pressure 

Bell  &  Howell 

1000-02 

9302 

■  ;-v- 

Transducer 

Digital 

Hewlett  Packard 

3466A 

1716A-18944 

Voltmeter 

DC  Power  Supply 

Hewlett  Packard 

6205B 

1140A14646 

Rotometer 

Brooks 

ShoRate 

USAF9175-83 

Travelling 

Gaertner 

slide 

2776P 

Microscope 

scope 

1211K 
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The  theore^cal  predictions  of  the  dynamic  response 
of  circular  crossAsection  pneumatic  transmission  lines  were 
compared  with  experimental  results.  Tests  were  conducted 
in  blocked  line  situations  as  well  as  in  both  laminar  and 
turbulent  mean  flow  conditions.  The  lines  used  in  the  study 
were  24  in.  in  length  with  inside  diameters  of  0.041,  0.119, 
and  0.195  in.  Tests  were  conducted  over  the  frequency 
range  20fjl000  Hz  and  over  a  range  of  mean  flow  Reynolds 
numbers  from  0-^16000. 

Theoretical  solutions  were  obtained  from  computer 
simulations  using  a  laminar  flow  model  and  two  turbulent 
flow  models.  These  models  were  developed  from  the  Navier- 
Stokes  equation  and  include  components  to  account  for 
viscosity  and  radial  heat  transfer  in  the  fluid. 

An  examination  of  the  turbulent  models  was  conducted 
based  on  a  nondimensional  frequency  and  attentuation.  This 
led  to  modification  of  one  of  the  turbulent  models  to  include 
a  frequency  dependent  term  in  the  conductance. 

The  effect  of  flow  conditions  upstream  of  the  test 
line  and  the  effect  of  orifice  design  on  the  dynamic 
response  were  studied. 

Experimental  results  showed  good  agreement  with 
results  of  earlier  studies.  The  use  of  a  bypass  device 
upstream  of  the  test  line  was  found  to  have  no  effect  on 
line  response. 

The  reconfiguration  of  the  orifice  devices  to  more 
closely  match  theoretical  models  produced  no  noticeable 
change  in  the  results, 
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